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PREFACE 

This final report as issued by the Applications Division of the 
Environmental Research Institute of Michigan (ERIM) under National Aero- 
nautics and Space Administration (NASA) contract NAS3-22442 for the 
Lewis Research Center (LeRC) covers the contract period from April 1, 
1980 through February 28, 1981. The technical representative for the 
contract officer was Mr. Thom A. Coney of LeRC. The Principal Investi- 
gator was Fred J. Tanis with important contributions to the technical 
program made by David R. Lyzenga, Glenn Davis, and Robert Dye. This 
research was conducted by the Applications Division under the direction 
of Mr. Donald S. Lowe. 

This contract involves developing algorithms to map selected con- 
stituent concentrations in Great Lakes waters from the Coastal Zone 
Color Scanner (CZCS). The approach is based upon the inherent optical 
characteristics of Great Lakes waters. 

A number of institutions and universities are involved in the 
project and are organized as the Great Lakes Experimental Team (GLET). 
This report covers ERIM's activities in the project during Phase I of an 
anticipated two phase program. 
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1 

INTRODUCTION 

The Great Lakes Experimental Team (GLET) is conducting a series of 
experiments in the Great Lakes region designed to evaluate the applica- 
tion of the Nimbus G Coastal Zone Color Scanner (CZCS). Potential uses 
foreseen include assessment of trophic status, verification and spatial 
refinement of whole lake models, and observation of temporal and spatial 
dynamics of phytoplankton. Currently members of the NOAA Nimbus Experi- 
mental Team (NET) are developing chlorophyll and sediment algorithms 
largely to be applied, to the open ocean [1]. Preliminary examination 
of these algorithms indicates they have limited applicability to the 
Great Lakes. Concentrations and compositional differences of suspended 
materials along with atmospheric aerosol variants are expected to exhi- 
bit important differences from the marine environment and result in 
additional complexities. The focus of the present program is the devel- 
opment and testing of atmospheric and water algorithms appropriate to 
the Great Lakes as well as evaluation of existing algorithms developed 
for the marine environment. 

1.1 STATEMENT OF THE PROBLEM 

The quantification of substances in Great Lakes waters by satellite 
visible radiometry is dependent on a thorough understanding of the radi- 
ative transfer processes in the atmosphere, at the waters surface, and 
in the water column itself. It has been well established that the 
content of water, be it particulate or dissolved substances affects the 
apparent color. By sensing color with a high signal to noise ratio in 
narrow spectral bands CZCS provides a means of looking at the water 
content which has been heretofore unavailable from satellite data. 
Since the air and water effects are coupled to the CZCS radiometric 
data, removal of atmospheric effects becomes critical to the success of 
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the Great Lakes verification. Once removed the radiance which is scat- 
tered upward from beneath the surface can be observed clearly by the 
satellite. Effectively the radiance reaching the satellite from the 
lake surface amounts to only five percent of the total radiation and 
consequently ninety-five percent of the radiation received is from 
atmospheric backscatter and surface reflectance. Furthermore, the var- 
iation in radiance at the satellite due to change in constituent concen- 
tration are on the order of one percent while the variation due to 
atmospheric changes can be considerably higher. The spatially varying 
atmospheric component is due principally to aerosol scattering. The 
significance of the atmospheric problem for a water target has been 
demonstrated by Hovis and Lung [2] and more recently by Quenzel and 
Kaestner [3] who compared the variability of the atmosphere with the 
reflected light from phytoplankton suspensions. Thus unless atmospheric 
effects can be negated resolving quantitative information on water con- 
stituents it is considered to be most difficult. Baring elimination of 
the atmospheric effects the water problem requires understanding how the 
inherent optical properties relate to the measured quantities of chloro- 
phyll-a pigment, phytoplankton cell count, suspended solids, and 
dissolved organics. Previous algorithms relied heavily on the availa- 
bility of extensive surface truth [4,5,6]. For this study algorithms 
are sought which can be based on optical properties specific to the 
Great Lakes and which reduce the present requirements for extensive 
surface truth. 

1.2 PROJECT GOALS AND ERIM TASKS 

In order to be acceptable to the Great Lakes user community CZCS 
algorithms must, in our estimation, meet at least two general criteria. 
First, the algorithms must be able to predict accurately surface concen- 
trations of chlorophyl 1-a pigments and suspended sediment over widely 
varying ranges and do so with little or no ancillary measurement data. 
Second, they must be capable of making predictions over water masses 
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which exhibit spatial variation in atmospheric haze and surface concen- 
tration. 

The ERIM participation in the Great Lakes experiments involves two 
phases. The first covers the period from April 1980 through February 
1981 and is the subject of this report. This first phase has involved 
development of computer software to process CZCS tapes received from the 
NASA Goddard Space Flight Center (GSFC), collection of surface truth 
measurement data in connection with CZCS overflights of the Great Lakes, 
and formulation of preliminary algorithms. The second phase will 
involve development and testing of specific atmospheric and water compu- 
ter algorithms for CZCS. The water algorithms will be based upon radia- 
tive transfer theory and measured optical properties of Great Lakes 
waters. Existing atmospheric models will be tested, including some 
recently developed models, using surface truth and low altitude aircraft 
measurements made during the 1980 GLET experiments. These models in 
turn will be utilized in the development of operational algorithms for 
removing atmospheric effects from CZCS data without the direct use of a 
large number of in situ measurements of atmospheric optical depth. Our 
approach is to remove spatially variable components using properties of 
the data itself. It is anticipated that both the algorithm development 
by Gordon [7] and that by a group at the Scripps Institute of Oceanogra- 
phy [8] will be tested for their suitability to the Great Lakes atmos- 
pheric environment. 

1.3 PROJECT BACKGROUND 

The scientific objective of CZCS is to determine water constituents 
quantitatively and to carry out such measurements over large areas which 
are not possible or practical to be obtained with surface ship investi- 
gations. Currently the Nimbus experimental team (CZCS-NET) is investi- 
gating CZCS capabilities to quantify material suspended or dissolved in 
the water. These validation studies are concentrating on the ocean 


3 



APPLICATIONS DIVISION 



environment. The present study is similar in design to the NET investi- 
gations but the focus is on a freshwater environment. A number of 
institutions, research centers, and universities plan to participate 
various aspects of the program. In addition to LeRC, which has led the 
current effort, participants in the Great Lakes Experiment include ERIM, 
Canada Center for Inland Waters, NOAA Great Lakes Environmental Research 
Laboratory, ERA Grosse He, University of Minnesota, University of Wis- 
consin, and others. These participants have a wide variety of back- 
grounds and capabilities which can be applied to the project. While it 
is understood that LeRC will not be able to participate fully in subse- 
quent program phases it is anticipated that they will maintain an active 
interest in the project and promote the continuity of the GLET. 

1.4 SUMMER 1980 GREAT LAKES EXPERIMENTS 

During 1980 a number of surface truth measurements were made coin- 
cident with CZCS overflights of the Great Lakes. Experiments were con- 
ducted at three principal locations in the Great Lakes: western Lake 

Erie, Duluth area of western Lake Superior, and the Grand Haven area in 
Lake Michigan. All of these experiments were designed to gather neces- 
sary validation data and optical properties specific to the Great Lakes 
waters. Measurements made included the following: 

(1) Aircraft flights were made by NASA LeRC F-106 aircraft fitted 
with the Ocean Color Scanner at altitudes of 500 and 41,000 
feet. 

(2) Water samples were gathered by the University of Michigan GLRD 
and the NOAA GLERL and subsequently analyzed for chlorophyll 
pigments and suspended solids. 

(3) NASA LeRC made various surface ship radiometric measurements. 
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(4) NASA Langley Research Center deployed a mobile optical labora- 
tory to Cleveland. Optical parameters including absorption, 
beam attenuation, and scattering were measured on selected 
water samples. 

(5) The Naval Oceans System Center conducted in-situ submersible 
radiometer measurements of subsurface downwelling and upwel- 
ling irradiance at selected sites in Lake Erie. 

A total of twenty separate sites were sampled in the Great Lakes, 
fourteen of which were made in connection with CZCS overflights. Only a 
very small portion of the above measurements have been analyzed to date 
but all analyses are expected to be completed in phase II of the 
program. 
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2 

DEVELOPMENT OF A GEOMETRIC CORRECTION 
ALGORITHM FOR CZCS 

The objective of this task was the modification of existing soft- 
ware and the development of a new scanner model which together will 
permit transformation to CZCS line and pixel coordinates into earth 
latitude and longitude. CZCS scanning geometry including variable tilt 
angle for Great Lakes viewing was combined with ground control points to 
derive the appropriate transformation matrix. Landsat geometric proces- 
sing programs were adapted and modified as necessary to process avail- 
able CZCS taped data for the Great Lakes. The initial accuracy goal was 
set to one pixel in each direction. In this task our efforts included 
investigation of geometric correction based on the 77 tie points per 
line, development of a CZCS scanner model, generation of mapping poly- 
nomials, modification of resampling software, and processing of an 
example image. 

The geometric correction of an image results from two operations. 
First, mapping polynomials are generated that define the transformation 
from the raw original image to the corrected image. Second, the correc- 
ted image is created from the uncorrected image using these mapping pol- 
ynomials. Two fifth-order, twenty-one term polynomials were used in 
this process, one for each dimension of the image. These polynomials 
define the transformation which makes the corrected image conform to a 
given map projection as well as adjust for viewing distortions such as 
satellite position, satellite motion, and earth motion. Two approaches 
were considered in the present study; a geometric regression analysis 
approach, and an orbit modelling approach. Both of these approaches are 
based on extensive experience with Landsat image correction. Thus the 
basic techniques were extended in the present effort to accommodate the 
Coastal Zone Color Scanner. 


PRECEDING PAGE BLANK NOT FILMED 
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2.1 INVESTIGATION OF REGRESSION MODEL TECHNIQUES 

Regression analysis was undertaken as the first approach to geome- 
tric correction. While this approach is considered to be straightfor- 
ward it has the disadvantage of requiring extensive ground control 
points which makes it time consuming. In order to get satisfactory 
results, fifty or more points should be taken for each scene. The first 
image considered was taken on November 8, 1978. An image analyst selec- 
ted forty-six ground control points of which forty-one were found to be 
suitable. A geometric regression analysis of these points produced 
mapping polynomials which could predict the location of these points 
with standard errors of 517 meters in the horizontal dimension and 553 
meters in the vertical direction. Since the pixel size of the CZCS is 
nominally 825 by 825 meters, derived mapping polynomials are estimated 
to be accurate to within one pixel. These results are comparable to 
those typically obtained with Landsat processing. 

Selecting ground control points is a lengthy process so an alterna- 
tive method was sought to correct the scene. An attempt was made to use 
the ephemeris data that accompanies each image tape. Specifically, 
anchor points are included that describe the geographic position at 77 
locations on each scan line. A two hundred point sample distributed 
throughout the scene was selected for testing. Geometric regression of 
these points showed errors greater than 13,000 meters in both dimen- 
sions. A similar set of points selected from a second scene taken on 
May 8, 1979 produced errors of the same magnitude. Subsequent to this 
experiment we learned that later versions of the GSFC processing algori- 
thm had improved the accuracy of the anchor points considerably. For- 
tunately, we were able to obtain a copy of the May 8 scene with the 
improved anchor point values. However, this image was the only one 
available under the new version [9] and thus it was used exclusively for 
purposes of testing the geometric correction algorithms. A geometric 
regression was performed on a 429 point sample of anchor points which 
was in turn used to produce a set of mapping polynomials. Using the 
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mapping polynomials derived from the anchor points the ground control 
points were predicted with RMS errors of 4068 meters in the horizontal 
dimension and 4425 meters in the vertical dimension. This five pixel 
accuracy is within that claimed by GSFC for the anchor point reference 
system. Geometric regression of a set of sixty-two ground control 
points taken from the same image showed errors of 557 meters in the 
horizontal and 893 meters in the vertical. 

An attempt was made to obtain additional scenes with which to veri- 
fy these results. Two scenes were obtained from the Lewis Research 
Center for this purpose. The first, April 17, 1979 had extensive cloud 
cover and it was impossible to obtain an adequate number of ground 
control points. The second, June 20, 1979 was centered in an area east 
and south of Lake Erie. Most of the image covered the Atlantic Ocean 
and again ground control points in the Great Lakes area were insuffici- 
ent. 


2.2 COASTAL ZONE COLOR SCANNER MODEL 

The scanner model developed for CZCS is based upon our experience 
with Landsat and in its present form takes each image control point in 
turn and projects it to the earth's surface. The line number for each 
point is used to interpolate for the latitude, longitude, and altitude 
from the values supplied with the tape reference data for each scene. 
The point number is used to calculate the mirror scan angle, which 
together with the reported tilt angle determines the scanner line-of- 
sight vector in spacecraft coordinates. 

A series of rotations through the angles of roll, pitch, heading, 
latitude, and longitude plus a translation provide the transformation to 
earth centered coordinates. The intersection of the line-of-sight 
vector with the surface of the eplisoid is then derived and converted to 
latitude and longitude. These coordinates are compared with the corres- 
ponding values obtained from map data, and the discrepancies minimized 
by successive refinement of the estimated roll, pitch, and heading. The 
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present model permits manual intervention with the operator who can 
supply needed refinements to the latitude and longitude parameters. It 
is anticipated that future model versions will eliminate the need for 
operator interaction and model adjustments. 

Once unsatisfactory control points have been eliminated the refined 
attitude and location data are combined with the model. The cartogra- 
phic projection is then selected and subjected to polynomial regression 
analysis which yields in turn the coefficients to a pair of twenty-one 
term, fifth degree polynomials. These polynomials provide an approxi- 
mating transformation from cartographic coordinates to the original 
image coordinates. 

The uses to which a scanner model of the imaging system for the 
CZCS are twofold. First, the image control points and their correspond- 
ing map control points may be easily evaluated for consistency with 
other points and any outliers rejected. Second, the coefficients for 
global mapping polynomials used in the resampling process can be derived 
by a fit to the model rather than to the points themselves. This 
feature permits the use of a much smaller number of points than would be 
needed for simple regression. 

2.3 GENERATION OF MAPPING POLYNOMIALS 

An arbitrary set of points were selected from the test image and 
used together with the scanner model to derive a set of mapping poly- 
nomials. A number of scanner variables including satellite attitude and 
position will influence the model results. Values of latitude and 
longitude are converted to line and pixel location in the resulting 
image via the equations that define the desired map projection. 

The process used to resample a corrected image works in the reverse 
direction. For each pixel location in the resulting image the program 
calculates the corresponding location in the original image. This 
process implies reversing the projection equations in order to derive 
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the desired latitude and longitude information. However, operation of 
the scanner model in reverse so as to select original pixel locations 
which correspond to given location in the corrected image is most diffi- 
cult and unwieldly. Alternatively, the pair of twenty-one term poly- 
nomials is generated to satisfy this mapping requirement. One polynomi- 
al describes the east-west position and the other the north-south 
location. The selected arbitrary image points mentioned above are used 
to generate coefficients for each term in the polynomial. The complete 
set of coefficients defines the mapping from the original image to the 
correct projected image. 

Mapping coefficients are generated by step-wise regression in which 
the correlation matrix is calculated relating each term to its ability 
to predict the location of the point in the original image. A regres- 
sion coefficient is calculated for the most influential term. This term 
is in turn removed from the matrix and the step regression continued 
until all terms that contribute predictive capability are included in 
the coefficient matrix. By utilizing points selected throughout the 
image file derived mapping coefficients are applicable over the entire 
pixel range in the corrected image. RMS errors are also calculated by 
the polynomial generation software for predicted locations in the 
original image. In the resampling process the correct image line and 
pixel number is translated to a location in the original image. The 
nearest neighbor pixel is then copied to the corrected image file. 

2.4 RESAMPLING OF THE CZCS MAY 8, 1979 IMAGERY 

Using the scanner model software fifty-four ground control points 
of the May 8, 1979 Great Lakes image file were tested for consistency 
with model parameters. Results indicated north-south RMS errors of 693 
meters and east-west errors of 1338 meters respectively. Based upon 
existing satellite attitude information the scanner model was used to 
generate latitude and longitude positions from which the appropriate 
polynomial coefficients could be derived. The polynomial prediction 
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errors were found to be 55 meters in the north-south direction and 583 
meters in the east-west direction. Combining these results leads to an 
expected total error of 748 meters in the north-south and 1921 meters in 
the east-west or approximately 1.0 and 2.5 pixels, respectively in the 
original image. Figure 1 shows the original and resampled images for 
the Great Lakes portion of the CZCS data file. This area includes all 
of the lakes except Lake Superior and the upper most portions of Lakes 
Huron and Michigan which are included in the next CZCS data frame. 

Unfortunately portions of Lakes Michigan and Huron are obscured by 
cloud cover. Lakes Erie and Ontario are essential cloud free with the 
exception of a thin covering over the western basin of Lake Erie. The 
resampled image was made using a polyconic projection and an arbitrary 
pixel size of 500 by 500 meters. So as to verify the accuracy of the 
corrected image twenty additional ground control points were selected 
and compared to those predicted by the projection polynomials. Table 1 
shows the results of this comparison. The maximum difference occurred 
for the sixth test point which was found to be 3000 meters in the north- 
south direction and 7500 meters in the east-west direction. Test pixels 
located near the center of the original image and principal meridian 
showed, on the other hand, minimal errors. For example the second test 
pixel had a north- south error of 500 meters and an east-west error of 
2500 meters. The mean error was estimated to be 200 meters in the 
north-south direction and 1875 meters in the east-west direction. A 
listing of CZCS geometric correction programs developed for the 
PDP-11/70 computer facility are contained in Appendix A. 
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Original CZCS Image 

Figure la. Original and Polyconic Projected CZCS Thermal Band Images of the 
Great Lakes for May 8, 1979. 
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Polyconic Projected Image 

Figure 1b. Original and Polyconic Projected CZCS Thermal Band Images of the 
Great Lakes for May 8, 1979. 




TABLE 1 

GROUND CONTROL LOCATIONS 

Predicted Actual Difference 


Location 

Row 

Col umn 

Row 

Col umn 

Row 

Col umn 

Toledo, Ohio 

1028 

1274 

1029 

1278 

1 

4 

Detroit, Michigan 

808 

1182 

806 

1191 

-2 

9 

Detroit, Michigan 

843 

1122 

841 

1131 

-2 

9 

Flint, Michigan 

692 

1381 

690 

1390 

-2 

9 

Flint, Michigan 

688 

1163 

685 

1174 

-3 

9 

Tawas City, Michigan 

311 

1150 

305 

1165 

-6 

15 

Cleveland, Ohio 

1039 

1564 

1039 

1567 

0 

3 

Erie, Pennsylvania 

801 

1495 

801 

1500 

0 

5 

Buffalo, New York 

731 

1806 

731 

1808 

0 

2 

Buffalo, New York 

865 

1827 

865 

1829 

0 

2 

Toronto, Ontario 

646 

1788 

647 

1793 

1 

5 

Toronto, Ontario 

658 

1895 

659 

1897 

1 

2 

Toronto, Ontario 

563 

1804 

563 

1808 

0 

4 

Elmira, New York 

785 

2204 

786 

2203 

1 

-1 

Elmira, New York 

715 

2308 

715 

2307 

0 

-1 

Rochester, New York 

635 

2129 

636 

2129 

1 

0 

Rochester, New York 

462 

2290 

463 

2290 

1 

0 

Kingston, Ontario 

438 

2103 

439 

2103 

1 

0 

Kingston, Ontario 

365 

2321 

365 

2321 

0 

0 

Kingston, Ontario 

419 

2222 

419 

2221 

0 

-1 


Mean Difference -0.41 3.75 

Standard Deviation 1.79 4.44 
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3 

USE OF OPTICAL MODELS TO DEVELOP CZCS 
CHLOROPHYLL AND SUSPENDED SEDIMENT ALGORITHMS 

A requirement fundamental to the validation of CZCS for the Great 
Lakes is the development of a working algorithm which can transform the 
satellite measured radiances into surface concentrations of chlorophyll 
and suspended sediment. While the ERIM task defined for current study 
involves development of water algorithms, these algorithms cannot in our 
estimation be attempted without some examination of and experimentation 
with atmospheric components. Thus, while we were able to place emphasis 
on certain water aspects of algorithm development, our approach has con- 
sidered radiative transfer in the atmosphere. Our efforts to date have 
involved extensive use of statistical and model simulation techniques. 
Attempts to test candidate algorithms on real CZCS data have been lim- 
ited because very few scenes of the Great Lakes were available and no 
CZCS tapes which correspond to the 1980 summer experiments are expected 
until mid 1981. 

3.1 APPLICABILITY OF EXISTING ALGORITHMS 

The removal of atmospheric effects is a necessary prerequisite for 
all remote sensing applications. Atmospheric effects are especially 
important in CZCS data for the following reasons. 

1. The large swath width of the CZCS implies a large atmospheric 
variability due to simple considerations of scale as well as 
scan angle variations. 

2. The inherent radiance of the water is low causing path radiance 
effects to be relatively more important than over land. 

3. The CZCS includes wavelengths in the blue region of the spec- 
trum where atmospheric effects predominate. 
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The Study of atmospheric effects in CZCS data can be broken down 
into two levels. First, one can attempt to develop and/or validate 
radiative transfer models by making careful measurements of the relevant 
atmospheric parameters and of the radiance at the surface, and comparing 
the radiance measured by the satellite with that calculated from the 
model. Studies of this kind have been carried out for aircraft data at 
various altitudes [10] and for CZCS data over the Gulf of Mexico [11]. 
Aircraft studies have resulted in fairly good agreement between model 
predictions and measurements, although there is a discrepancy at large 
angles which is thought to be due to surface reflected skylight [10]. 
Previous studies with CZCS data have encountered some difficulty in 
obtaining agreement between model predictions and measurements [11]. 
One possible explanation for this difficulty is the effect of scattering 
from adjacent land areas which some studies have indicated to be of the 
same order of magnitude as the directly scattered path radiance [12]. 
In its studies ERIM will test existing models with aircraft measurements 
made during last summer's GLET experiments and with data obtained in the 
Gulf of Mexico Experiment [11]. 

A second kind of atmospheric study involves the development of 
operational algorithms for removing atmospheric effects from CZCS data 
without the use of an unreasonable number of in situ atmospheric measur- 
ements. The primary goal of these studies is to remove the variable 
component of the atmospheric effect using some property of the data 
itself to obtain the necessary correction parameters. One such algori- 
thm was developed by Gordon [7] and applied to CZCS data over the Gulf 
of Mexico [10]. However, the formulation of Gordon's algorithm involves 
the assumption of zero intrinsic radiance in the 679 mm band which is 
not met in some parts of the Great Lakes and other coastal areas. Modi- 
fications to this algorithm have been made by a group at Scripps Insti- 
tution of Oceanography for conditions occurring in Pacific coastal 
waters, but there is some doubt that their assumptions would hold in the 
Great Lakes. In addition to testing these algorithms, new directions 
have been pursued in the development of more suitable approaches. 
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3.2 OPTICAL MODELS AND RADIATIVE TRANSFER THEORY 

Most existing algorithms are based upon empirical relationships 
between constituent concentrations and remotely sensed radiances, and 
are generally valid for the limited range of environmental conditions 
under which these relationships were derived. In order to systematical- 
ly approach the evaluation of existing algorithms or the development of 
new algorithms, it is necessary to understand the relationships among 
the observed quantities on a more fundamental level. The study of these 
relationships is conveniently divided into two phases. The first phase 
deals with the inherent optical properties of the water aand atmospheric 
constituents, and the second phase deals with the large-scale radiative 
transfer processes which relate these inherent optical properties to the 
radiances measured by the satellite. 

3.2.1 OPTICAL MODELS 

A full description of the optical properties of a passive (non- 
emitting) medium include the absorption coefficient and the volume scat- 
tering function. For most remote sensing purposes, however, it is not 
necessary to specify the complete volume scattering function. Commonly 
two parameters describing the scattering properties are considered: the 
total scattering coefficient (which is the integral of the volume scat- 
tering function over all angles), and the back-scattering coefficient 
(which is the integral of the volume scattering function over the 
angular range of 90° to 180° from the incident direction). 

It is generally assumed that the inherent optical properties are 
linear functions of the concentrations of the various constituents of 
the medium. Thus, for natural bodies of water we can write 
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N 

Bb=Bb„+fBbC, (1) 

i=1 - 
N . 

b = M b^C. 
i=l 

where a , Bb , and b are the absorption, backscattering, and total 

WWW /V /\ /V 

scattering coefficients of pure water; a^ , Bb^. , and b^ are the 
absorption, backscattering, and total scattering cross sections (i.e., 
the coefficients for unit concentration) of constituent i, and is the 
concentration of constituent i. These are clearly approximations, since, 
the actual optical properties may depend upon factors other than the 
concentrations. The scattering properties of suspended particulates, 
for example, depends upon the size distribution as well as the mass per 
unit volume. Some preliminary studies of the effect of size distribu- 
tion are described in section 3.3. 

Accepting the linear model for the optical properties described 
above, there are several possible approaches to the determination of the. 
absorption and scattering cross sections for each constituent. One 
approach would be to attempt to isolate each constituent and measure the 
optical properties for different concentrations of that constituent. 
Isolation can obviously be accomplished more easily with some types of 
constituents than with others. The danger of this approach is that by 
artificially changing the sample, the measured conditions might not be 
representative of naturally occurring waters. 

A second approach for determining the optical cross sections is to 
measure the optical properties and the constituent concentrations Of a 
large number of diverse natural samples, and performing a statistical 
analysis (i.e., a multiple linear regression) of the results. This 
approach avoids the possibility of encountering unrealistic conditions, 
but is subject to biases introduced by natural correlations between 
various constituents, for example between phytoplankton and sediments 
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associated with phytoplankton nutrients. Examples of this approach 
include the work of Bukata et al [13] and the preliminary analysis 
presented in section 3.3 of this report. 

It should be pointed out that some of the optical properties, i.e. , 
the scattering cross sections, can be calculated from electromagnetic 
theory if the size distribution, shapes, and indices of refraction of 
the particles are known. Although it is instructive to make such calcu- 
lations in order to determine the sensitivity of the cross sections to 
the size distribution, for example, there is no particular advantage to 
the exclusive use of this method since it requires measurements which 
are at least as difficult as the direct measurement of the optical 
properties. 

3.2.2 RADIATIVE TRANSFER THEORY 

The radiance at any point in the ocean/atmospheric system is pro- 
vided, in principle, by the solution of the radiative transfer equation 
with the proper boundary conditions and the proper optical properties 
specified at each point. Unfortunately, an exact solution of this equa- 
tion is not obtainable in closed form for even the simplest geometry. 
Two general types of approach are, therefore, possible. The first 
approach is to develop approximate solutions which give reasonably ac- 
curate results over the range of conditions encountered. The advantage 
of this approach is that results can be obtained at a relatively low 
cost, and that considerable insight can be gained by merely examining 
the functional form of the solutions. 

The alternative approach is to develop exact numerical solutions of 
the radiative transfer equation. Several types of numerical solutions 
exist, perhaps the most prominent for this application being the Matrix 
Operator and Monte Carlo methods. The Monte Carlo method is especially 
powerful because of its ability to incorporate complex boundary condi- 
tions (e.g., at the water surface) and spatially variable optical pro- 
perties. The advantage of the numerical approach is, of course, the 
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accuracy of the results. The disadvantages are the cost of obtaining 
these results and the fact that insight can be gained into the nature of 
the solutions only by examining a large number of cases with different 
input parameters. 

The preliminary modeling work done on this project has involved a 
combination of the approaches described above. For the purpose of form- 
ulating candidate algorithms it has been useful to consider a simplified 
model in which the water radiance is calculated using Gordon's Power 
Series Approximation, and the atmospheric effects are described by a 
model which assumes constant transmittance and a path radiance which is 
a linear function of the aerosol optical depth. According to this 
model, the radiance at the satellite is given by 

L = Lp(x,) . TEp„()() (2) 


where 

= path radiance (assumed to be a linear function of the 

P o 

aerosol optical depth t ) 

a 

T = atmospheric transmittance (assumed constant) 

E = irradiance at water surface (assumed constant) 


P„(x) 

X 


^ ^wl ^w2 ^ 

Zirn^ (y + y^) 

= (c.f. previous definitions of a and b in text) 


^wl’^w2 


= water surface transmittance for incoming and outgoing 
light 


Pq, y = cosine of solar zenith angle and observation angle, 
respectively 


n = index of refraction of water 
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The algorithms developed using this model are described in sections 
3.4 of this report. For the purpose of evaluating these and other 
algorithms, a more comprehensive simulation model was also developed. 
This model incorporates the Monte Carlo calculations of the subsurface 
water reflectance using the following power series expansion [14] 


Pw(x) = — ■ 2 ^ [.0001 + .3244x + .1425x^ + .1308x^] 

7T n 

where x is defined as previously. Atmospheric effects are calculated 
using the QSS model for the path radiance, and the double-delta model 
[11] for the irradiance and sky radiance. The surface-reflected sky 
radiance is included, assuming a nominal surface reflectance of 2.0 
percent. The atmospheric state is described in terms of the horizontal 
visibility using the relationships developed by Elterman [15] to calcu- 
late the optical depth. The results of this simulation model are 
presented in section 3.5 of this report. 

3.3 PRELIMINARY ANALYSIS OF OPTICAL PROPERTIES FOR THE GREAT LAKES 

During late July 1980 twenty-one individual samples gathered by 
LeRC were analyzed by the NASA/Langley research Center (LaRC) portable 
laboratory stationed at the LeRC flight facility. Samples were flown in 
on the same day as collection and usually received by the LaRC staff 
within four hours for analysis. Underwater optical properties were mea- 
sured in vitro with identical spectral range (400-800 nm) and intervals 
(50 nm). These properties included absorption, beam attenuation, and 
volume scattering. Three separate instruments were used to measure 
these parameters, (1) a combination beam attenuation and small angle 
scattering meter (SASM, 0 = 0.379, 0.751, 1.49°) developed by LaRC and 
patterned after the Scripps Institution of Oceanography ALSCAT instru- 
ment; (2) a Brice Phoenix (BP) scattering meter modified to accommodate 
large angle measurements (25° ^ 155°); and (3) the LaRC spectral ab- 

sorption coefficient instrument (SPACI) [16]. Standard errors for these 
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instruments are reported to be as follows: (1) for the SASM less than 
5% a, and less than 12% e(9); and (3) for the SPACI less than 10% a. 

The optical measurements made by LaRC during the 1980 summer exper- 
iments were for lake samples. In this case the measured optical proper- 
ties pertain to the particular mix of constituents in the lake sample. 
If sufficient number of measurements are made in this manner and if the 
principal constituents present are known then multiple regression tech- 
niques can be used to derive the optical cross sections for a common 
constituent. While the present data are considered limited for this 
purpose a preliminary set of optical properties were derived for the 
Great Lakes. 

Of the twenty-one optical data sets taken, three Lake Erie samples 
contained sufficient quantities of sediment to saturate the optical mea- 
surement instruments. Samples collected from Green Bay in Lake Michigan 
and from Western Lake Superior were found to have distinctive local 
optical properties. Four of the six samples collected from the Grand 
Haven area were essentially sediment-free and the presence of very low 
concentrations of phytoplankton made absorption measurements difficult. 
Attempts to include these samples in the regression analysis have so far 
not been productive. 

The best regression results were obtained when nine samples from 
Lake Erie were combined with two samples from Lake Michigan. Several 
regression models were formulated and tested against the above selected 
optical measurement sets. These models were based upon the available 
surface truth sampling data which included Secchi depth, surface 
temperature, chlorophyll-a pigment, phaeophytins, and residue (total, 
ashed, and volatile) [17]. Of these parameters chlorophyll-a, chloro- 
phyll-a plus phaeophytins, total residue and ashed residue were selected 
for regression. Each of the candidate models contained two components 
and a constant which includes absorption or scattering for pure water. 
Models involving ashed residue produced generally better statistics than 
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those with total residue. Each regression model consists of four equa- 
tions pertaining to the backscatter cross section. In each case the 
four equations correspond to four CZCS wavelengths (443, 520, 550, 670 
nm). The optical model considered for these analyses describes the 
surface water mass to be a combination of pure water (w) , unique 
organics as represented by chlorophyl 1-a (chi) and phaeophytin (pp) 
concentration and unique inorganics as represented by the measurement of 
suspended minerals (sm). The two component model equations are written 
as 


a(^) = ^ constant a(x) 

Bb(x) = Bb^(x) + xBb(s^l^(x) + yBbgj^(x) + constant Bb(x) 

where x and y are the concentrations of chlorophyll a and suspended 
minerals (ashed weight) respectively, and a and Bb are the absorption 
and backscatter cross sections. 

Optical cross sections as derived from these regression analyses 
are given in Table 2. These optical cross sections are the only such 
data available, to our knowledge for the Great Lakes. Note that the a^ 
or Bb are included in the constant term given for each analysis. As 
shown in Table 2, two preliminary optical models were derived based upon 
chlorophyll-a and chlorophyll + phaeophytin concentrations, respective- 
ly. Also shown for comparison are values of four and five component 
Lake Ontario models [18] which were derived indirectly from apparent 
rather than inherent optical measurements. The four component model 
included chlorophyll, suspended sediment, pure water, and dissolved 
organics. The five component model has an additional term for non 
living organics which includes detritus. The dissolved organics term 

accounts for the presence of yellow substance which was found in the 
Ontario study to be about 2 mg/£ and fairly constant throughout the 
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TABLE 2 

OPTICAL CROSS SECTIONS FOR GREAT LAKES WATER QUALITY MODELS 


(1) Great Lakes 1980 Preliminary Models 


Model la 


Model lb 


Wavelength (nm) 

2 


Constants 

443 

.01620 

.0764 

.3088 

520 

.00836 

.0636 

.2484 

550 

.00529 

.0577 

.3034 

670 

.00450 

.0556 

.3897 


BbcHL(>"^/mg) 


443 

.000152 

.0312 

.0424 

520 

.000372 

.0284 

.0370 

550 

.000469 

.0287 

.0232 

670 

.000428 

.0250 

.0197 





443 

.01142 

.07290 

.3794 

520 

.00599 

.06178 

. .2813 

550 

.00384 

.05658 

.3220 

670 

.00323 

.05468 

.4072 


®'>CHL.pp<"'^/"’9) 



443 

.000043 

.03117 

.04590 

520 

.000213 

.02832 

.0408 

550 

.000273 

.02856 

.0277 

670 

.000254 

.02490 

.0237 


(2) Lake Ontario Five Component Model 


acHL('"'^/"’9) 




Constant 


443 

.0354 

.0557 

.020 

520 

.0240 

.0281 

.028 

550 

.0173 

.0185 

.037 

670 

.0100 

.0225 

.370 


Multiple 

Regression 

Coefficient 

.961 

.952 

.945 

.942 


.900 

.911 

.911 

.913 


.970 

.963 

.954 

.947 


.890 

.910 

.910 

.912 
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2 2 
/mg) Constant 


443 

.00199 

.0328 

0 

520 

.00182 

.0474 

0 

550 

.00241 

.0525 

0 

670 

.00175 

.0333 

0 


(3) Lake Ontario Four 

Component 

Model 


2 

a__(m^/g) 
sm' ^ 

Constant 

443 

.0343 

.0557 

.185 

520 

.0232 

.0281 

.119 

550 

.0173 

.0185 

.122 

670 

.0105 

.0225 

.388 


2 

BbcHiCni /n»g) 


Constant 

443 

.00163 

.0328 

.0010 

520 

.00153 

.0474 

-.0058 

550 

.00202 

.0525 

-.0099 

670 

.00156 

.0333 

-.0080 
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lake. Since for the present study no effort was made to analyze the 
samples for a non-living organics or dissolved organics term they were 
not included in our preliminary model except as part of the constant 
term. 


The absorption cross sections as derived from the 1980 optical mea- 
surement sets are for chlorophyll, about half of those derived for the 
Lake Ontario models. Derived chlorophyll backscatter coefficients were 
found to be only one tenth of those obtained by the Lake Ontario study. 
Chlorophyll absorption and scattering cross sections as obtained by 
regression are considerably less than those reported elsewhere in the 
literature [17]. On the other hand, derived optical cross sections for 
the suspended sediment are comparable to the Lake Ontario values. 

All of the above models have utilized phaeophytin free chlorophyll 
determinations. Since chlorophyll and phaeophytins have similar absorp- 
tion properties they cannot be readily distinguished by CZCS. There- 
fore, it seemed appropriate to combine these determinations in the above 
regression analyses and determine cross sections for chlorophyll plus 
phaeophytin. The resulting cross sections as given in Table 2 model lb 
are smaller for the chlorophyll term and larger for the constant term 
than in the phaeophytin free analysis of model la. Since the direction 
of these changes are opposite to our expectations, additional investiga- 
tion into the optical properties of phaeophytins is warranted. 

In addition to the statistical analyses of the collected samples 
described above a brief investigation was made of an indirect technique 
for measuring scattering and absorption properties of suspended sediment 
in Lake Erie. The approach involved making a series of physical and 
optical measurements on a single turbid water sample at various time 
intervals. A portion of the suspended sediment was allowed to settle 
out between each measurement. Optical parameters measured by the LaRC 
portable laboratory included the beam attenuation coefficient and the 
volume scattering function at 90°. Both of these measurements were made 
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at a single wavelength of 550 nm. Physical measurements included the 
particle size distribution (using a Coultor counter) and the total 
suspended solids concentration. There were two objectives to these 
measurements. First the change in suspended solids concentration could 
be related to the corresponding change in the optical properties from 
which optical cross sections could be obtained. Second, Mie particle 
scattering theory could be used to calculate the same optical properties 
based upon the change in particle size distribution. Unfortunately, the 
total suspended solids measurements were found to be unreliable, either 
because of sampling or measurement errors, and were not included in the 
analysis. The optical measurements are summarized in Table 3, and the 
particle size measurements are shown in Figure 2 (panel a). It should 
be noted that Coulter counter measurements are difficult and subject to 
considerable error in the particle size range necessary for studies of 
this kind. 

Because of the unreliability of the total suspended solids measure- 
ments, it was not possible to derive the optical cross sections directly 
from this set of measurements. Instead, the data set was used to test 
the feasibility of the approach of calculating the optical properties 
with Mie scattering theory. Of course, these calculations also requires 
a knowledge of the index of refraction of the particles. Since no 
measurements of the index of refraction were made, an average value of 
1.1 was assumed. The extinction efficiency was calculated from the Mie 
theory using J.V. Dave's algorithm [19] for particle diameters between 
1.0 and 25.4 ym. The portion of the total extinction coefficient 
arising from this particle size range was then obtained for each sample 
by multiplying by the size distribution function and integrating 
(summing) over this range. For particle diameters less than 1 micron, 
the extinction efficiency was calculated using Van de Hulst's 
approximation [20]. 

^ext = 2 - ^ sin p + -^ (1-cos p) (6) 

^ P 


\ 
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TABLE 3 

OPTICAL MEASUREMENTS MADE DURING SETTLING EXPERIMENT 


Sample 

Settl i ng 
Time (hrs) 

a(550 nm) 

e(90°, 550 nm) 

1 

0 

23.7 m'^ 

0.1380 m"^ sr'^ 

2 

1.5 

19.3 m“^ 

— 

3 

19.0 

8.7 m‘^ 

0.0685 m"^ sr"^ 

4 

23.5 

7.0 m"^ 

0.0617 m“^ sr"^ 

5 

40.0 

4.5 m"^ 

0.0319 m'^ sr"^ 
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where 

4TT(n-l)r 
P ' A 

(n is the index of refraction and x is the wavelength). The contribu- 
tion to the total extinction coefficient from particles less than 1 
micron in diameter is than given by 


ext 




(7) 


For this particle size range, the size distribution was assumed to have 
the form 


dn 

dr 



( 8 ) 


where c is chosen to fit the measurement of r = 0.5 ym. The extinction 
coefficients calculated from this procedure are shown in Table 4, and 
are compared with the measured extinction coefficients in Figure 2 
(panel b). In view of all the uncertainties and assumption required to 
make these calculations the agreement shown is surprisingly good. While 
these results are obviously insufficient to validate the approach they 
do suggest that further controlled studies of this type would be benefi- 
cial. If a technique could be derived for universal Coulter Counter 
measurements instead of the present elaborate and not readily available 
optical measurements the overall benefits of the present program would 
be substantial. 


3.4 DISCUSSION OF ATMOSPHERIC EFFECTS AND CORRECTIONS 


A necessary preliminary step before attempting to extract any in- 
formation about the water itself is to remove the effects of atmospheric 
variations from the measured radiances. Atmospheric effects are 
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TABLE 4 

EXTINCTION COEFFICIENTS CALCULATED FROM MIE SCATTERING THEORY 


Sample 

Contribution From 
r < 0.5 ym 

Contribution From 
r > 0.5 ym 

Total 

ext^' 

1 

7.2 

18.3 

25.5 

2 

4.8 

14.1 

18.9 

3 

2.5 

6.7 

9.2 

4 

1.0 

4.3 

5.3 

5 

0.7 

3.2 

3.9 
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particularly important in CZCS data because of the large swath width and 
the band placement. The development of radiative transfer models for 
the atmosphere is an important step in understanding this problem, but 
ultimately the goal must be to remove the atmospheric effects on a 
point-by-point basis with the aid of only a very limited number of 
external measurements. 

The first attempt to formulate such an atmospheric correction 
algorithm for CZCS data was made by Gordon [7]. This algorithm is based 
upon the assumption that the total radiance in a suitable wavelength 
band may be interpreted as an index of the atmospheric state (i.e., the 
aerosol optical depth) and thus used to correct the atmospheric effects 
in the other bands. In this formulation it is assumed that the water 
radiance in the 670 nm band is zero or negligible compared to the path 
radiance, an assumption which is valid in clear ocean waters but is 
frequently violated in more turbid coastal waters, including the Great 
Lakes. In fact, our modeling study has shown the 670 nm band displays 
the greatest sensitivity to sediment concentration. Therefore, although 
the algorithm gives apparently good results in ocean areas [1] it cannot 
be applied directly to the Great Lakes. 

The assumptions about the atmosphere in Gordon's algorithm seem to 
be valid for a reasonably wide range of atmospheric variations. These 
assumptions are essentially the same as those listed for the simplified 
model described in section 3.2, namely: (1) the path radiances in the 
various wavelength bands are linearly related to each other, and (2) the 
atmospheric transmittance changes relatively slowly and may be assumed 
constant. Under these assumptions one can define a large class of 
linear combinations. 


^i 


N 

I 

j=l 




(9) 
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which are independent of atmospheric variations. Gordon's algorithm is 
one member of this class, but is not necessarily the optimum one. The 
condition for the X^'s to be independent of atmospheric variations can 
be written as 


N 

T A.. 
^ 1J 
j=l 


9L. 


= 


N 

I 

j=l 


A.. A. 
1J J 


= 0 


( 10 ) 


where i=l...N-l, L„. is the path radiance in band j and t is the aero- 

PJ o 

sol optical depth. The coefficients A., may be interpreted as compo- 

->■ ' J 

nents of the vector A^ , and the above condition viewed as the require- 
ment that these vectors be perpendicular to the vector 


(A-| , ^2 > 



( 11 ) 


If we require in addition that the vectors A^... A^j_j be orthonormal, 
i .e. , 


A. 



0, i j 

1 , i = j 


( 12 ) 


then the transformation (9) has the properties of a projection onto a 
hyperplane perpendicular to A|^, and the projected variables are linearly 
independent of each other. Although this would seem to be a desirable 
condition, the actual benefits of this procedure, as opposed to Gordon's 
procedure for example (which does not satisfy equation (12), have not 
yet been demonstrated. Work has begun on the evaluation of this 
procedure using the simulation model described in section 3.2, but a 
full comparison with Gordon's algorithm has not yet been completed. 

It should be noted that the direction of the vector Aj^ can be 
determined empirically from the CZCS data itself, if an area of variable 
haze over a uniform water background can be located in the image. The 
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direction cosines are obtained readily by a principal components 
analysis of the radiances observed over such an area. 

As discussed previously in section 3.2 the signal variants due 
to spatial changes in the aerosol content (haze) are linear. As 
discussed by Gordon in his correction method the aerosol contribution at 
one wavelength is approximately proportional to the other wavelengths. 
In CZCS four channel signal space, atmospheric variation is visualized as 
a vector oriented by these proportions and offset from the origin by the 
presence of water constituents and atmospheric transmittance effects. By 
comparison the water variants represent separate orientations for each 
constituent and generally much smaller than the atmospheric haze vector. 
Thus for open clear waters of the Great Lakes the observed variations are 
due essentially to atmospheric haze and system noise. 

In order to explore the haze phenomenon several segments from two 
available CZCS images were examined. Each segment (200-1000 pixels) was 
first scaled to radiance and then analyzed for principal components. 
Two available CZCS scenes were selected for analysis: Great Lakes May 

8, 1979 and Gulf of Mexico, November 9, 1978. Results of these analyses 
are shown in Table 5. Also shown is the principal component for 
atmospheric variation as derived from simulations using the preliminary 
atmospheric model discussed in section 3.2. Atmospheric vectors derived 
by this analysis show strong similarity in orientation. Differences in 
orientation are likely due to cloud effects and variations in water 
constituents. Radiometric changes due to scan angle effect may also 
account for some of the observed differences. Excellent agreement 
(1.1°) was found between the combination of Lake Huron data sets and the 
theoretical orientation of the haze vector. In general, the orientation 
of the haze vector for the Gulf of Mexico data sets were similar to 
those for Lake Huron but with less consistency and lower accounting of 
the percent of total variance. Clouds and variation of chlorophyll and 
suspended materials could possibly account for the observed variability 
in principal components. 
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TABLE 5 

FIRST PRINCIPAL COMPONENTS DERIVED FROM 
CZCS SATELLITE MEASURED RADIANCES 


Sample Location 

First Principal 
Component 

Percent of 

Total 

Variance 

Georgian Bay, Lake Huron (1) 

(.563, .539, .543, .313) 

94.8% 

Georgian Bay, Lake Huron (2) 

(.440, .499, .565, .488) 

94.2 

Southern Lake Huron 

(.435, .577, .594, .355) 

95.4 

Combination Set, Lake Huron 

(.426, .5520, .547, .500) 

95.5 

Gulf of Mexico (2) 

(.596, .580, .467, .298) 

98.3 

Gulf of Mexico (3) 

(.487, .548, .538, .416) 

94.8 

Gulf of Mexico (4) 

(.559, .597, .465, .339) 

91.7 

Gulf of Mexico (5) 

(.604, .679, .330, .256) 

85.1 

Combination Set, Gulf of Mexico 

(.521, .579, .518, .353) 

87.4 

Atmospheric Model 

(.426, .510, .524, .533) 

100.0 


CZCS Imagery 



Lake Huron: May 8, 1979 Orbit 2715 

Gulf of Mexico: November 9, 1978 Orbit 227 
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An analysis was made to observe if the haze vector was present and 
significant in apparent low haze area. Several samples were selected in 
the May 8 image from a test area south of Nova Scotia in the open ocean 
and at least 160 kilometers from the US mainland. The image appeared to 
be free of haze. Nevertheless the haze vector appeared as the principal 
component in each sample and subsample and with as few as fifteen 
pixels. These analyses further demonstrate the need to account for the 
atmospheric variants even under the clearest of conditions. 

3.5 PRELIMINARY CZCS WATER ALGORITHMS 

The atmospheric, interface, and subsurface water reflectance models 
discussed in section 3.2 were used to calculate expected satellite 
radiances for a variety of water masses at each of the CZCS wavelengths. 
The primary input for these calculations were the optical cross section 
data as described per the three optical models in section 3.3. Differ- 
ent water masses were simulated by varying the concentrations of 
chlorophyll and suspended mineral concentrations. In the case of the 
Lake Ontario five component model the level of non-living detrital 
material was taken as 2.0 mg/£ and dissolved organics at 2.5 mg/£. 
These levels are similar to those measured in the Lake Ontario study 
[13]. Presently we have no measurement data to support the 
representativeness of these values to Lakes Erie or Michigan. 

Having made these assumptions the subsurface irradiance reflectance 
can be readily calculated for each CZCS band (443, 520, 550, 670 nm) as 
a function of the concentration of chlorophyll and suspended minerals. 
The spectral characteristics of the irradiance reflectance function can 
be depicted with iso-concentration curves for each pair of wavelengths. 
Figures 3, 4, and 5 show calculated subsurface reflectance for each of 
the two Lake Ontario models and the preliminary 1980 optical model, 
respectively. Each of the four panels of each figure has nine curves of 
increasing suspended mineral concentration and constant value of 
chlorophyll pigment concentration (0.0, 1.0, 2.0, 5.0 10.0, 20.0, 50.0, 
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Figure 3. Subsurface reflectance (percent) at CZCS wavelengths (443 nm, 

520 nm, 550 nm, 670 nm) as predicted by the Lake Ontario 5-component 
model [18]. Each figure has nine parametric curves of increasing 
suspended mineral concentration (0.0-50.0 mg/2.) with constant 
values of chlorophyll pigment concentration (0.0, 1.0, 2.0, 5.0, 

10.0, 20.0, 50.0, 100, 200 yg/£). Each figure also contains four para- 
metric curves of increasing chlorophyll a (0.0-200.0yg/£) at constant 
values of suspended mineral concentration (0.0, 1.0, 10.0, 

50.0 mg/£) . 
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Figure 4. Subsurface reflectance (percent) at CZCS wavelengths (443 nm, 

p20 nm, 5b0 nm, 670 nm) as predicted by the Lake Ontario 4- component 
model [18]. Each figure has nine parametric curves of increasing 
suspended mineral concentration (0.0-50.0 mg/£) with constant 
values of chlorophyll pigment concentration (0.0, 1.0, 2.0, 5.0, 

10.0, 20,0, 50.0, 100, 200 pg/£). Each figure also contains fourpara- 
metric curves of increasing chlorophyll a (0.0-200.0 yg/£) at constant 
values of suspended mineral concentration (0.0, 1.0, 10.0, 

50.0 mg/£). 
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Figure 5. Subsurface reflectance (percent) at CZCS wavelengths (443 nm, 

520 nm, 550 nm, 670 nm) as predicted by the Preliminary 1980 
Optical Model. Each figure has nine parametric curves of increasing 
suspended mineral concentration (0.0-50.0 mg/£) with constant 
values of chlorophyll pigment concentration (0.0, 1.0, 2.0, 5.0, 

10.0, 20.0, 50.0, 100, 200 yg/il). Each figure also contains four para- 
metric curves of increasing chlorophyll a (0.0-200.0 yg/£) at constant 
values of suspended mineral concentration (0.0, 1.0, 10.0, 

50.0 mg/£) . 
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100, and 200 pg/i). These curves tend to converge to a point beyond the 
calculated range. One could refer to this ideal point as the "point of 
all sediment". Each figure also contains four curves of increasing 
chlorophyll with constant values of sediment (0.0, 1.0, 10.0, and 50.0 
mg/£). These latter iso-concentration curves tend to converge in 
several of the panels to a point on the 670 nm axis. This ideal could 
be referred to as the "point of all chlorophyll". Each panel of these 
figures is a projection and together suggest the reflectance space is a 
three dimensional hyperplane and nearly perpendicular to the 
R(550)/R(520) plane. 

For the five component Lake Ontario model the constant optical 
cross sections are very small relative to those for chlorophyll and sed- 
iment. As a result the sensitivity of reflectance to changes in concen- 
tration is large. By comparison the Lake Ontario four component model 
has slightly smaller cross sections and a larger constant term. Conse- 
quently the panel figures are slightly smaller. The effect of the 
relatively small optical cross sections of the preliminary 1980 optical 
model with a large constant term is shown by the small magnitudes of 
change in each panel of Figure 5. The difficulties of using this latter 
optical model to predict concentrations is apparent. The above figures also 
show generally that reflectance is more sensitive to changes in sediment 
(as mg/i.) than chlorophyll (as yg/Jl) which is consistent with their 
relative optical cross sections. 

As discussed in the previous section, the spectral changes observed 
in CZCS data due to chlorophyll and sediment will be influenced by the 
presence of atmospheric variants. A principal component analysis of 
pure chlorophyll and sediment data indicated spectral orientation with 
angular separations from the pure atmospheric vector of 26.1° and 18.7°, 
respectively. Thus it seems apparent that the atmospheric variants are 
indeed coupled to those we wish to determine in the water. Unless there 
is some way to separate the atmospheric and water components by spatial 
filtering it seems appropriate to remove this influence by projecting 
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these reflectance data along the atmospheric vector. This projection 
reduces the four dimensional reflectance space to one of three 
dimensions which is free of any atmospheric influence. Figure 6 shows 
the projected space for each of the three optical models under consid- 
eration. The panel figures on the right are nearly a planer view of the 
depicted leaf like projected structure. Thus corresponding projected 
axis hold promise for decifering the water components. Since there are 
an infinite set of transformations which will project the four dimen- 
sional reflectance space into three dimensions the panel figures shown 
can be rotated to any orientation. This feature may provide a means to 
later optimize candidate algorithms. 

Thus far our efforts to obtain chlorophyll and sediment algorithms 
have utilized the above projection technique. Using non-linear 
regression techniques algorithm prediction equations were derived as 
third order, second degree polynomials with nine terms. These equations 
have the following form: 

f(conc) = Cj^(Rj^) + C 2 (R 2 ) ^ 3 (^ 3 ) ^ ^ 

II I I I ( 13 ) 

+ C7 Ri + Cg R2 + Cg Rg + CjQ 


where the R's refer to the three axis of projection and the C's are the 
multiple regression coefficients. 

Using the model described in section 3.2 and the optical properties 
of Table 2 data sets of simulated CZCS radiances were generated for 
water quality conditions similar to those that exist in the Great Lakes. 
Several simulation data sets were generated for each model. These 
included each of the following types: 

(1) Variable chlorophyll, fixed sediment, fixed haze, zero system 
noise; 
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(2) Fixed chlorophyll, variable sediment, fixed haze, zero system 
noise; 

(3) Variable chlorophyll, variable sediment, variable haze 
condition, standard system noise level. 

The simulations included the effects; of system noise by adding a 
normally distributed random variable with standard deviation equal to 
the mean radiance divided by the signal -to-noise ratio. The signal -to- 
noise ratio used were those reported by Gordon for CZCS [1]. The pre- 
liminary algorithms derived from these simulatipn data sets use as a 
first step the atmospheric projection procedure described in section 
3.4. The orientation of the atmospheric vector was that given in Table 
5. Data set (1) as described above was then used to derive a prediction 
function of the form given in equation (13). Similarly data set (2) was 
used to derive a corresponding equation for sediment. The third type of 
data sets (3) were used subsequently to test performance of the derived 
prediction equations. 

Results of these analyses are summarized as follows: 

1. Sediment equations for the Lake Ontario models were able to 
predict sediment concentrations under conditions of variable 
haze and system noise to within approximately 50% over wide 
ranges of concentrations (1 to 20 mg/Ji) and to within 30% over 
narrow ranges (1 to 5 mg/£). 

2. The sediment equation for the preliminary 1980 optical model 
were able to predict sediment concentration to less than 50% 
only after the signal -to-noise ratio was doubled. 

3. None of the regression equations for chlorophyll were capable 
of predicting chlorophyll to within 50% under standard system 
noise conditions. When the noise was reduced by four times the 
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Lake Ontario optical model predictions fell within the 50% 
error range. The presence of sediment, as might be expected, 
was found to have a deteriorating effect on the prediction 
algorithms for chlorophyll. 

4. Satisfactory results were obtained with the preliminary 1980 
optical model only when the presence of sediment was reduced to 
very small quantities (« 1 mg/t) and the system noise was 
virtually eliminated. 

While we are encouraged by these simulation results additional 
investigations of this type will be needed under Phase II in order to 
produce satisfactory algorithms. The simulation techniques developed 
in the present study will become a powerful analytical tool for 
investigating and evaluating the applicability of Phase II CZCS 
algorithms to the Great Lakes. 

Attempts to apply derived chlorophyll and sediment algorithms to 
Great Lakes CZCS data sets were not viable because of the lack of 
surface truth. However, some preliminary analysis was performed on the 
May 8, 1979 scene. The derived atmospheric vector compared well with 
that obtained from the atmospheric model as discussed in the previous 
section. 

The processing portions of the final algorithm is anticipated to be 
able to classify the image and separate it into a land and water file. 
The water file, which will contain the surface area of the Great Lakes, 
will be then radiometrically calibrated in the first four bands. In 
turn, the water file will be transformed to one suitable for constituent 
determinations by projecting the data along the atmospheric vector. 
Applying the sediment and chlorophyll prediction equations to the 
projected file will produce the desired pixel by pixel concentrations 
maps. These concentration data files would then be recombined in the 
computer with the land mass file to produce a final map product. 


46 



APPLICATIONS DIVISION 


2erjm 


A first attempt at performing some of these manipulations was made 
for the May 8, 1979 scene. While complete processing was beyond the 
scope of the present program phase we were able to separate the image 
into land and water files, color slice the water file, and recontiine the 
separated image as shown in Figure 7. 
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SUMMARY AND CONCLUSIONS 

Much of the work accomplished to date is preliminary to the valida- 
tion of CZCS in the Great Lakes. During the first phase of an anticipa- 
ted two phase program, efforts were initiated and directed toward 
development of the necessary algorithms and supporting processing 
software which will allow the transformation of CZCS images of the Great 
Lakes into maps depicting concentrations of chlorophyll -a and sediment. 
The second phase will involve analyses of the CZCS imagery and surface 
truth measurements collected during the 1980 surmier experiments. 

Preliminary examination of existing CZCS atmospheric correction 
algorithms developed by NOAA for the open ocean indicates that the 
assumptions required are not valid for much of the Great Lakes area. 
Work has begun on the development of new atmospheric correction algori- 
thms which are appropriate to the Great Lakes. While these algorithms 
appear to be promising, they have not as yet been thoroughly tested with 
actual CZCS data. 

A preliminary optical model was derived from the LaRC optical mea- 
surements made for Great Lakes waters. Derived chlorophyll -a cross 
sections were found to be less than those reported with the Lake Ontario 
models [18] and elsewhere in the literature. 

Efforts to derive chlorophyll and suspended sediment algorithms 
were based upon simulations of the preliminary optical model and the 
Lake Ontario models. Simulations included system noise and atmospheric 
variants as calculated using existing models. Multiple non-linear 
regression techniques were used to derive water quality prediction 
equations from Great Lakes CZCS simulation data. While results produced 
are encouraging they are thus far incomplete because of the lack of 
sufficient optical data and appropriate CZCS images. Suspended sediment 
algorithms were found to be able to predict sediment concentrations with 
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single pixel accuracy to within 50% of the true value over the range 
(l-20mg/A) for all optical models under consideration. Chlorophyll on 
the other hand was found to be more difficult to predict because of its 
smaller optical cross section. A 50% prediction accuracy could only be 
obtained after substantial reduction was made to the system signal to 
noise ratios. Furthermore this improvement was only realized with the 
Lake Ontario optical models. Satisfactory chlorophyll predictions could 
not be obtained using the preliminary optical model for the Great Lakes. 
This result was anticipated in part since the derived chlorophyll 
optical cross sections were much smaller than 
expected. 

A second activity of the current work involved development of a 
geometric correction algorithm for CZCS. A scanner model specific to 
CZCS was developed which accounts for image distorting scanner and 
satellite motions. This model was used in turn to generate mapping 
polynomials that define the transformation from the original image to 
one configured in a polyconic projection. 

Actually two approaches were investigated in the present study to 
obtain these mapping polynomials; geometric regression and orbit model- 
ing. Based on a single available CZCS scene for the Great Lakes a geo- 
metric regression of anchor points produced mapping polynomials which 
predicted the location of ground control points with RMS errors of 
approximately five pixels in both the horizontal and vertical direc- 
tions. By comparison the scanner model produced RMS errors of less than 
one pixel in the horizontal and 1.5 pixels in the vertical directions. 
Thus the scanner model approach is presently considered to be superior 
to exclusive use of image anchor points. 

While some minor modifications to the scanner model are anticipated 
as additional imagery is acquired the software package to provide CZCS 
geometric correction is essentially complete. 
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APPENDIX A 

CZCS GEOMETRIC CORRECTION SOFTWARE 

This appendix contains FORTRAN IV listings of four programs as 
developed for the PDP-11/70 computer system under the present contract. 
These programs form the basis of ERIM's present capability to transform 
NASA/GSFC CZCS image files into an image with desired metric qualities. 
The four programs include: (1) the CZCS Scanner Model, (2) CZCS Image 

Ground Control Program, (3) the Mapping Projection Polynomial Generation 
Program, and (4) an adapted Nearest Neighbor Resampling Program. These 
programs do not constitute a stand alone capability but instead require 
selected supporting software from ERIM's Earth Resources Data Center 
(ERDC) operating system. Available documentation and operating 
instruction can be obtained by request from ERDC. 
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.330 PEAL*> (iwOHK 


■ < 7 .A ' r - . 

3-;-r lOGiCAL*! LHF( 1 P73C1) ,HFN(1«} ,ST CIPT ,0PIJE (1 1S0S5 

3 9^1 LOC JC al * 1 S I'M 1 PI , HT ( 5P) , COORDS ( a) , SWAP , LOAD .77 7 . 

aci? 7 BYT E ' ANCH'IP , SE C f ai 

a u’ “ ■ equivalence t YE ar ,lBE (1 ni , f SECNDS , LHF f ? 1 ) ) ,7' 0 A Y , LRF 1 1 9)5 

aP0 ' ■ EQUJ valence CT U.T,LBF n 7) ), CTILTID.LPF (19) ) 

a >00 I VAt ENLE ( OEUE .OPI.lF? ) , ( I MGS 7 , L8F ( ?5 1 ) , C GMT SE C , SE C ( 1 ) ) 

ELiiiTVALEMCF (LBFfn.PE.CNBR) 

as0 "" EDuTV'AlEnCE 7 f LATLOM,C(iOROS) , ( RESULT,C0nR0Sl7 77 ' 7/ 

aA07 7 77 ^7 7 77 7 _ __ . 7. ._7. 77.77 77777 

'a 70 ~’~'-c ' 7_ "y “ y 7 7 " .7..77r 

a'“0" c 

aopi ■ r y'fVPEN E I LEs,’ SET ' CONST ants 7 _. . . 7. 7 _ . 77. 77 77 77 7 7 

^00 C " " ' 

A10 ASS IGN~ ' TO I99aPR 

s?0 t;n ID Roaon 7 __ _ 7_ 1_ ' _ _ 7 7 _ 

' ss'd ’’”''99a'0? ’cONi iNut ' ■ 

'sau ■ r y 7 'y 

330 C' GET E IPST' nadir 'values 

3A0 'C ' * 
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APPLICATIONS DIVISION 


SFO C7CSMKr.,FTN ?h-FEH-8L PAfiE 2 

! ! I t t t 1. 1 

'?7i7i ASsSIGN 908("? TO I99etH(!| 

SP.Ti GO TO 99Ane« 

59PI 9~96>? 'CONT iNUE 3 ' T 7 ' ■ ■ 

C -- ^ 1’ 

[ i p " c — — j— - J-- Z I ' Z’Z 

c 

I" ' r r7ET FI«'ST GBdi,JP OF' ATTITUDE' values 

f>o'{5 C‘ ’ " ■ 

['■ ' ' t.s?l ' _C" ASSJW' 99J5i?''^^^ I9930B IZ~ I IT I 7 '7 

FE? C" GO TO 99500 


! 

t.7 0 

C9950 

? CONTIMU: 





PE0 

c 

- 




L...... 

E90 


THGHEG a IMGEDG 





71? (71 

c 





! 

710 

C 






7P0 

9ti0 

CALL RkAOfUNIT.LHF, 12780, NAP, 

IC) 



1 

7 5? 


IE" ( NAR.NE . I ?78Pi) ' GO TO 950 





7U0 

C' "■ 

IF-"t .NOT “ft OAD.OR. ANCHOR) ) GO 




I 

750 


TO 900 




7e0 

C'" "" 





[■■■ 

77? 


"SWAP i' LPFM T 




r" 

7P0- 


"IRFfn 9 LFF(2) ' " "■ 




7 90 


LFF(21 a SWAP '■ 





f 00 

"C 






P 1 0 

"0 

"type 910,HfcCNUW,RECNSR 





P 

,0910 

FOPBAT ( I 5,08) 





® ? ' 

"C ' ' " 








' BECNPW a 971 - ( RECNBR /16 ) 

♦ 1 





*’>'1 


I 

ADO ONE TO hake value 

ONE RFLATIVE 


a *.« 


J 

SUBTRACT FROM 971 TO 

invert 

image 


fl 7f? 

A. A'.’! " 

0 "" . 

_J 

DIVIDE BY 1A> TO SHIFT 

RIGHT 

a BITS 



r 

TYRE 910,PErM8R,»F.CNRR 







TILTIO a MtlOfT ll.no, 25S) ~ ■ 




■ 

9)0 


IF 1 TILTIO.NE.2) GO TO 66«l 





9?0 SNAP * LPFf) 7 ) 

(■"■ 950 "'LhF(l 7 ) » lBF()flr' 

940 , ■LR‘^ 1 1 P V »■ SWAP 

r" '9SO ' c ■■ ■ '■ 

' 9 f ,0 TF 'ITU.T.FO.TILTSV) GO TO 66 (? - 

f 970 T]|;TSV b tilt ■ 

9F0, C niv int' BY ?5S ■ TO USE LEFT half' OF “WORO WHICH IS' THE 'WHOl.E PART ‘ 

i 090 "0 OF THE NiiMREP 

‘■"_'j 000 RwnKY"'="FL(JAT( ' TILT ■/ "256 ') 

[■ I'H fli '"" “ TLTVAL'' HWOPK"' *' -,''3h7 ♦ '29',«7' ■ ■ ■ 

10?, 0 ' ■ 0 WPIif r5,*,H?V PWOWKjTLTVAL.RtCNflR "" 

^"1070 D'M?~"FORHaT f ' TILT s ',F 5 . 0 ,* ' TILT "fOEGT" « 'VF&Vfl",'' AT RECORD" 15 ) " 

f C < 1 0 - - g 0 - C'n N T T N I ) E - • 

r "■ ) 0 S 0 ""'c' 

' lApo c r " 

I "' 1070- TF r'VNOT,'ANrHnR)'"GO: TO 'psa " 

■ ' 1 c 

r '1090 c ■■ "PFT'LIME'OF anchor POINTS 

" 1 1 00 C ASSIGN 9090 ? TO 199900 “ - ' - 

<110 C GO ■ TO ■ 99000 

" 1 1 ?0 ' "'COOB/? CONTINUE ' : - 



rz 
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SEO 

1 1 5PI 
1 1 

r_'l 1 5 pi' 

■ 1 t ^ P! ■ 

r 11 7 PI 

1 1 HPI 

[’"’l 1 ‘3pi' 

" ! ? PI (i 

r~ 1^1 P! 

■ 

L ’ 1 ?ipi ' 

[""''1P5PI' 

i-jfrtr 

I " f ?7P| ' 

1 ? p 

; 

1 <BIP ■ 

L ' n 1 ' 
- 

[" 1 TJPI 

' II'JPI ■ 

f " l^’S'P 
V '» P ?l 

t 

[' 1 pi^n” 

ippipi" 

I"'" it) 1 P! ' 
\iiS-i' 

1 aaiii ■ 

f " 1 «spi" 
1 ahpi 

C.' i 

1 ap-ii 

i 1 « << PI 

1 

I " ! S ) '7 

IPI’P! 

r ' 1 S P PI ' 

isapi' 
isppi ■ 
iS^n' 

f ' i PTCI' 

1 sa'ji 

Q ' IS'-JP " 
“ 1 h CHI 

r" 1 1> 1 ii ■ 

1 

i "ibiei"' 
1 ii>ap 
1 ^S"l 

" 1 f.hl7 ' 

i“ l(S7il' 

Vspc* ' 


Pei-FEH-81 l2ia«lP7 
>1.; I i 


PAGE 


C7CSMRG.FTW 
1 1 

C 

c 

C ^ 

bea~'TF T,NiiT,LnAbT''Gb to'^9 
II l._ bn' hRt) BA140 ■ r ' NOVI “ ' 

I„IST««T « RfeRIHG ♦ (BAND - "n ' * NE 

fc90 CALL MOVE (L_3F(sf4Wn , 1 fOPUFfBANOI ,NOVL,ME,03 

c " 

b ’ I. II w » T T F ( K B ,'9 p en ■' ( 1. H F ( I ) , I •> e G I M G , fl 6 6 ) , ( 0 B U F { J 3 , J ■ 1 , 

3f 6 ( lx,033 ,3X) 3 3 


.30), 


. c * L4_ . J. 7 P C N B 

I'bn ^IPl^Pn ' II^ J ' " 

950 ' COn't I'lijUfi ' 


IF_(L0A(M caul close (1)_ 

TA pLfJSE (UNIT • 33 


c 

c " 

C"- ■■~GFT SECONin SET OF NlADiw' values 

‘assign 99fl03 TO I 99 P,! 3 PI 

'Go''Yb' 99a0ci ' 

99 9'ii’3 ' CON 1 1 NiJF 

r 

"c "i ’■ 

C 

'c 'GET'SECnNO GPbijp OF ATTITUOE'IvALUES 

'C 'I''_I 

ASSl'CN 99?«?_T0 ''isb'pppi _'I^ 

'c "go 'to ■roPP’Y'J ' ■ “ 

'C9'9?0?_CnNf INUE' '"'I' 

C . - 

C' . 

c _ Copy naoips points to nadir file 

"r ' II'I“I^ _i 

__ A'g si'o'N 99piY?"TiY _ _I I! 

I'go rn" 99 'n 0 pi_’ ’• ^ 

990 P> CONTINUE “ - - 

"c" ' -I'l I 'll J 

C ■ I ■ '"'I' 

"■"“"■"ff ■"( ,NOTluIb 4 lb 3 ~jGO"^^^ 2 " "I "" 

C " ■" _ 

'C' " - 

'r ■ ' " 

c ■LbAo'''HeAbE> file' ' 

c" 

ASSIGN 99hb?2fo'IT99600'''^2''I'2"^^^^^ I 
go’ TO 99SH0 ' 

■ 99 S 0 > CONtiNUE ' '22 21 ' tr _ IL _ r 

C ' ' ' 


c 

5i'0' 


CLOSE I UNIT ■ a 3 

CALL re “ TunI T")“ 


56 
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SFO 

1 >.9C1 

I" ITIC'I 
i 7 >ei 
[ "l 
1 7 Ufi 

' 7 - ’ 

' i 

r’"l 77P 
"1 7«C! 

[■ i 79 '.l 

I H '^71 

( 1 P 1 ? 

1 P ?91 

i' ' 

' ‘ 1 ««« 
r tpspi 

' t MPP 
1 1 A7!7 

■ 1 AA'<I 
t flqp' 

1 9 </ 0 ' 
1 9 U' 

1 9?5! 



“■■ ■l 9a?! 

1 990 

1 

■"■ 1 970 

1 9P0 

1 900 

?O7i0 

?0 ) 0 

?0?0 
' ? 0 ! 0 
? 0(10 

■ ?oso 

? 0 9 0 
3070 

f ?090 
'? t 00 

?iia 
" 3 1 ?o 

■ ? t AO 
?1 00 
?1 so 

■ ? 1 <1 0 

70 

? 1«0 
? ) 90 
" ??oo 
i ' '■ P? ! 0 

??«0 


C7CSMH5.FTN PP-FHB-91 12l«at{»7 PACE 0 

I I i ....1 t t 1 


C99901 GO' Tb' 1 ' , • • 

99B01 r.O'TO l998ObV(*J9Ao2,9980'5) 

c J'" - 

C9970T GO' Tn t997jjj,- - 

9970i 'G'O to 199700', ■'t<»‘'7M8, 99707, 997387 ' 71'_1 . . - 

c 

998(Dr'G0"T0 I9o»,0(il, (99803) " ' ' 

c " 

r9950t' 'Gti TO j 99500, (99?0a, 99503, '99504, 99505)7'- 7 

c T 

99«or r,0" TO 1994 00 , (99 4 03 ) ‘ 

c; 

r.99J0l"GO' to 199300,(99303) : 

C 

099P0r'G0_JO_' 199300, (99303) 3 — - — 

C 9 9 1 PM G 0 _T 0 1 9 91 0 0 , ( 9 9 1 0 3 , 9 9 1 0 3 , 99 10 4 ) 

9 9 1 ■■■■GO" T O'" T 99 B PI 0",_( 99003) 7 — — - ■ - 

ill ***** '**'*'•*'*•'* * *****'*'** *r* •*♦**'**'**'*■*'*«’*'* **'*’***’*'*' 

C ' ■"■GET “AMCHMJ’OIMTS — ^ 

'C 9'9'9O0"'CO'N'T't «UE' - - 

'r' ■■■" 


I 




C 30MGITUDES' MUST BP 'cqNVEOTEO 'PROM 0 TO" 360 SYSTEM' TO +- 180 OES, 

'r'"" ^ ^ ■_ 

c '"7— * 

C E4C'H GE0GM*PMIC‘C00BnTNATE" OCCUPIES FOUR BYTES'; 

c""" Pit 31 Reppfsemts the sigm 

C BT)S 50-33 REPRESENT THE WHOLE DEGREES PORTION' OF THE NUMPER. “ 

'c ■■""Pits ?i-o represent the sev'en-digit decimal fraction of degrees. 

r ■■ 

C "STNE' s' LPF' (' WMC4NC ) .AND, ''300 1 SIGN IS ALWAYS POSITIVE 

r '"SINE s IS iGN'n, sine ) 

C ’_ 'I I I ■ — - ^ - ^ - 

c" ■"* ■; ■'■■' • 

f _ 'W'h'CPNT"* "o' ■■ '" ----- — - — - - - -- - —■ 

'r 

C 0 0 ■ 9 8 3" W H C A N C ” s ' 9 E G A N c , E N 0 A N c , 4 

"t'~ ■' ■ ■ ■ ■■ - 

'C'"'" 

C DEGREE '"i 'LPF (' WHCANC )",AND,"'"177 * " 

C OFGRFE ■ * ISMFT(0EGREK,3) 

C WORK S LBF ("WHCANC ♦ T ) ' , A NO , ' " 3 00“ 

C ^"WORK"* ' WORK 7 "84 ; \ ' SHIFT " 6"_UITS T0'"R1G«T"""™3 

, degree: + woRK ^ ' " _ ' " _ " ■ _ '_ " ■ 

C"" FR4C1 *LPF ( WHCANC ♦ 1) ,ANO, "077 
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SEQ CZCSMWG.FTN P 6 -EE 8-81 tPl««IPr PAGE 5 

f < i > 1 . t I 1 , . . . 

r. EflAC? «L 8 E ( WHCANC ♦ 2 ) .AND, " 37 ? 

PPfsVi c 

r ■" 2 ? 7 flj ■ ' 2 'c ” FJ 7 A C 3 1 ' *'l p F ■ J ■■ W H c A N C J+ 3 ) "Va NO , ' » 3 7 7 ' ' ’ ' 

??«oi _ C __ _ "■■■■_ -■ 

j'y ppq./ C ‘'jNlir «' NHCANC' 3 ’ _ 

. ^ -Z - : 

L j>._ TVPK oil ,kHCAnC, aBFf nVlBWHCANCjENOri iLBFViri 

,.f.OJ>MAT( • POINT ',I 5 , 2 X,» OtCf », afih,lX) _ 1 

L ^ ♦ ,/jJ«X,^OCTl '.flfOfejlX) _) 

? 3 a 0 C _ _ ■ 

[_ C FNAC « FPACl * "rwoife ♦’ FRAC 2 * TWOB ♦ FRAC 3 

? «*.? c 

C_'” 7 Z «ZflOAT ■( FRAC ) '/ ( TW 032 5 ~ 'J " "1 Z 

?<p.VI ' r. ~ _ ' ' ” - - - - - 

L '> 39 yl -J GEO' « "" TLO’aT f DFGRt e ) '♦ RWORK 1 Z Z Z. ZZ Z 

?ap«! ' ' c ' ’ ■' 

rz ■■■' CZZIZLJ^ ‘geo" ■ «'"( 3 «»B. - GEO V z 

’ c’ z” ~ z _ "Z" 

r P'J AS* ■ n 'tvPe; '‘RipVfJf&^EF.FRACl ,FRACP, FRAC?, FRAC, RWORKZTWOai ' ' 

CS 30 FORMAff a f fB. IX) , r lei.FlS, IP.FIB,?! ' ' 

I z 2 usp ‘ c ' , ' " zz ' z. z - _ _ _ 

' P'JbiA D TyWR.?'?, '■' 6 Ert' ■ ■' 

■ f ">070 cr^'Zformat fj X,' Fi 2 V») ■ 'ZZZZ’Zr * 'ZZ'-”' 

?«AC 1 c ’ ' 

r ?aR?lZ C w'hCp'nT "**'wHCPNf i' _ 

"" PSP 7 I c ‘POtNTSfNHCPNf) VgEO Z 

i”' ?sirti - Z — _^' 2 Z" ■ z Z’ 

"" r.RM? _continue'Z”* _ _ ' 1 Z _ Z - _■ Z _ Z" 

[' ■'"‘'r 'ZZZZ ZZ ZZ _Z 1 ' _ 

?^H!A t ■ 


[_■ ?«iS!>l ‘'Z C z* WRTTE’r.VRECNHR) f PblNfsTKy, K"» ^ 

?ShB pQi'ilATT. 

[Z r 7 ZZ c ■ ' ~ "z 

?s«M z, c z.._ 


?e'/p r Z .Z.I'X’yZT^" 


czz_z 

?^.3'7I ■ r_ _ZZ^6 cFO0rE' T o ''ExfRACT‘‘NA‘r>IR' TN 

?N«7i z”f^ ~ _.zz ' z ' zzi-.,” zizz -Z z z z .Zf 

Zc Z ~Z!Z'Z ..'.ZZZ1ZZ1..ZZZ''ZZZZZZZZZ .zz'zyzzzz, z,_ 


■' ‘>NbP RRflBy CdATlNUe Z.ZI' 1 _ 1 ZZ..ZZZ Z 

! pbtb' " c r " Z Z Z,..Z__Z 

C_"_Z 1 ” 

r~ r. ' ' get ~GNT' INAGE'EOGE^CiMGEnG) AND' LENGTH OF IMAGE (SPAN J 

'Z'PTW'A c'”' 

r ? 71 v» C _ SwItcm' BYTES and WHRPS TO" MATCH DEC CONVENT IONS 

"■? 7 ?(iZZ'Z Z- Z ' ■ 

r“ 2 7 3 n> ' ZZ„LypoZ'*?iT' 7 >N'<Z*-.J!i - . 

?7 «r‘ZZ.C - — y - - 

i SNAP ■ I.HFCWORK) > 

‘ ? 7 hp "■ “ Lrf ("worh' i s lbf (Wor‘k "'4 3 'Y ’ ..Z.Z'ZZ,^ ...Z 

IZ ?7 7 ?l L 6 F( work' j» 3 ! 

? 7 «B C . _ 


I " ?7R0 ' ' ' swap .■ LSFfWPRK ♦ ) ) ZZ Z.Z Z Z 

■ ■ t-RF C WORK * t ) ■. -L ZV 
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ORIGINAL PAGE ES 
OF POOR QUALITY 

APPLICATIONS DIVISION 


sso 

Pepc/i 

? s « i? ■ 

? fl 7 PI 

??oni 
> Q 0 '.1 ’ 
1 'll 
? J?SI 

?^»SPi 

r a <» ; Cl 

'■ ■ ? q f< CV" 

C 

TtPS*" 
[■■' iPM p) ' 
" Ti 1 ?pi' 

I " T PI 3^1 
Ti/iu?i' 
I"' 5 ''IS pi ' 

Tilfr'?" 

i 7 'a"' 

ipii^pi ■ 

F JPI9PI' 




CZCSMHG.PTN 

).. I. 


SPi 


LBPf WOBK ♦ ? ) 
CONTINUE 


?A-FCB-Ht {2144107 

J i I 


PARE 


• SWAP 


I ... 


l._ 


c 

b" 
■ r 

“c' 

■'C"'' 

c 

c 

■:c7 
c ' 
c" 

C' ■ 


IMCeOG ■ SECNOS ^ 

SPAN • iMliSZ 

IMRfcDG' « 1MGE05 ♦ SPAN-'J 
TYPt ■ . , * TPREOG • . iMCebC, 


SPAN »» 8 PAN 


'VALOf G 1 VtM bN' TAPE FOR SPAN 'DOESN'T APPE AR TO' BE 'REL I AHLE 
SO A OEFAULT SIZE OF 128 SECONDS IS USED, -- - • 

"sPAN''» 't2s00o''' 


'~'R)(IHACT time of GMT hlNuTe MARK t "'MINMRK') ■ 

■_ j i N S C R A M H L E G M T' j / 1 ? 0 A V 

offset it' HEGILT ♦ OFFOAY 

GMTSEC ■ "a 

■■ SECCa) ' « LPP'r PPFSET ) - 

SFCtn ' • LPF 'I offset V 1") 

GMT01 2 '•' GMTSEC ■ 


_c 

c" 

"f" 
■'C " 
-C'- 


I 


3 1 ■ 

["'■.31 50 ' 
■■■ 31 Ufi' 
[[' 3! S51' 

■ A1 

[■' 2.3170 
3 1 « ' 71 ; 

1 3 1 9171 ' 

saiTHi 
[' ■■" 3 P 1 0 ' 
3?aO' 
3?30 
3 P«PI' 
saSiTi 
3?fec7i 
3a7j! 
AaSi:' 

.3 a 90' 
3 5?"a 
331 i? 

3 3 a'^ 

3 3 3IA 
33 «b 
33S7) 
33 ''b 


EYfLOOE DAYS HEFORE CURRENT OAT 

"WORK s hOO" fGMTOia, 1 2) _ 

CONVERT TO MTLHSECOMOS"' 

'work s' work •' 2 58007100"' 


■ UNSCR AMHl £ CP T M ILL I SECONDS 0P ‘ 1 / 1 2 ' OA Y ' 

"OFFSET '■ BFCILT'+ OFFSEC ' --j " 

'CMTsSEC' » 0'" 

SEC' ( 3 T"» LEIF'C OFFSET) “ 2 

"SEC f'a 'l "a "IfiP f OFFSET ♦ l ") 

■'SEC f " l r s' LBF ( OFFSET ♦ 2 " 


r 


L 


L... 


AOb MlLLISECONns OF THIS TWO HOUR SEGMENT" " ' " " 

" MINMRK ' s'GMTSEC' '4''W0RK 

TYPE «, • WORK '.WORK,' GMTSEC '.GMTSECy* MJNMRK 

'22'' 2 Choose MI NUT F_M ARK inoexes 

"" RTIMF 3 "float r MInMHX '/■ E0000 ) 1 , 

GTBASFs MINT GT RASE, RTIME ) 

T r P E * 7 ■ ' G T B A 5 E ' , G T M A S E 


GET NADIR "VALUES "" 
no 200' WHCMIN s 1 ,2 
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SEQ 

CZCSMHG.FTN ?6«PEB-8t 12144107 PAGE 

7 




.1 



NAOCNT ■ NA0CN7 ♦ 1 




STA9T ■ flEGlLT ♦ GMTSZ ♦ (WHCMIN « n * 

POSSZ 

i 3391^ 


11FK3ET ■ NLONSZ ♦ START 



C 



1 38 1 Jl 

C 

GFT LON value 


3u?n 


ASSIGN 99706 TO 199700 


1 3a3<’> 


GO TO 99700 


ja«!3 

99706 

continue 


1 3as^ 

C 





NAOLONTNAOCNTl « ELOATC RESULT! * OPMR 


r 3a7Pi 

C 



3«P0 


OFFSET A NLATSZ ♦ START 


r 3a90 

c 

GET LAT value 


SSPipi 


ASSIGN 99707 TO 199700 




GO TO 99700 


3S?a 

99707 

continue 


r 3S3!i 

C 



3SaCl 


NAOLAT (NAOCNT! ■ FLOATt RESULT 1 * OPMR 


r 3ssf 

c 



35bl?i 


UFKSfeT ■ NALTSZ ♦ START 


r 3S7? 

■ C ■ 

GET ALT VALUE 


3'>«0 


ASSIGN 99708 TO I99T00 


r 3S90 


GO TO 99700 


3h'<VA 

99708 

CONTINUE 


1 3bl0 

3*»P0 


NA0*LT (NAOCNT! ■ RESULT 


\ 3b 33 

C . 



3ba0 


NADTIM(NA0CNT! ■ ( MINMRK ♦ ( (WMCMIN - 

1) A 60000 )} 

r 3bS0 


♦ f 60000 


3bb!? 

C 



1 3b70 

0 

Type «, » nadir lat.lon.tim, alt '.naolat (naocnt! ► 

3be0 

0 

♦ NAOLONI NAOCNT! , NADT I M (NAOCNT) , N AD ALT (N AOCNT! 

L. 

c . 



3700 

?00 

CONTINUE 


i 3710 

c 



37?0 

c 



L 

c 



3700 


GO TO 99801 


! 37*10 

• r * * * * 



3 7*>0 

c 



! 3 7 7 0 

c 



3700 

c 

OPEN FILES, SET CONSTANTS 


r 3790- 

c 



3800 

c 



r 3»I0 

99400 

continue 


38 00 


PI ■ 3. 1415927 


f 38 30 


OPMR « i.E-6 • 180 / PI i DECREES PER 

MILLIRADIAN 

^fl«0 


TWO.?? » ?, •• ?? 


r 38S0 


TW016 ■ ? ** 16 


3860 


TWOH ■ ? «• 8 


f 3870 


NPHATT • 0 1 


3«80 


NAOCNT 10 \ 

. 

f 3890 


GT6ASE ■ 999999. 


3900 

C 



f 3910 

c 



3980 


K B s S 




r;: 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


APPLICATIONS DIVISION 


V 


SEO 

CZCSMRG.FTN 

1 ,.,.1,,, 

28-FEH-8I 

>i 

12IA4IP7 

PACE 


IDR*1 

C 


1 

DRIVE 1 

B95P 

WRITE a<e,?en 





?Sl ■ format' (' C 7 CS'bA"TA TAPK MERGE VER ji, 0 ', 7 ) 

P T‘> 7 iX ^ 

■■ ■ ■ _pWRf Tf ■ {K 8 ',>P)“' ‘ ■_ ’7 

r Tbp JJNip*Nb7bENSlTV? _ ' 

ofl(«9i "r " “ ■ ' 


8 


t. 



[ C THTS' STATe"EMT HEAOS 'A DIAGNOSTIC' CONTROL ' VARIABLE 

ap?9t C J ''AnCi' i > LDAOS hoTH 'IMAG t AND ANCHOR POINT 'FILES ''' ' '7 ' 

[■ R5)?ci C ANS ■ « 'l LOAPS ONLV IMAGE FILE 

“ '4PUH c " "ANS "• B lOAPS ONLV "ANCHOR POINT FILE 

(" ■ UPISP ■■ READ 'fkp/iiPl' UNIT, OEN, ANS ; 

(IMbi.i RB' FORMAT 'LPIIPI ■ 

1 4Pi7pi c ''77 _ "I ZZ 

aeiSvA CALL-'INTT (UNIT,0EN,'2) 

[ '■"' call ' HfcWTNI) fUNIT) 

■■'411 pp •""■'C" 

f" All P IF r 'ans.nE.P,'anD,AnS','NE,I "')~00 TO 72 • 

'~ '"AIPP LOaD'''«"'.tkuf .' 

r A13P' HR! Te"tKB,«iP1 "■ ■■ " 

4 UP 50 forma t " r ' SSCEnE 'TlTtE^t O' : 

f" A15P PEAO' fk'H, *,pV ST'" 

4tSt» Bp ' FORMAT (i JOAl'y ' — 

7 'a! 70 "c ' 

«UP ' "WRrTE "'('kB'.'BSV"' 

I A 1 9P " FOR MAT" f ' SOH T VE"NUMBEPT'”'* I" 

"A3PP " RE'kl) (KB'.apy'TOR . 

“ A31P 'C - ^ ■■ r — 

"" A?3p WRITE "'(KH,'7P) ■ 

r A ? J P - 7 p F 0 R M A T' f » $ H E A 0 E » ' T I T L E t ' " ” , 

a?up PHAir 'f'Fh/fcp) 7" ” 

f '"'A?Sp- 't"'' :'■■ 

a ?bdi 7 ? COf<iTTNUE^ 7 

r A?7p c"'"’ ■ 

A?8P C ■"anchor"'* .FALSE'," 

l'''A?9’^ C IF (AnS. N fc.P. AND. NE.2) GO TO 7A 

'■ A3PP C 'A'NCHOR' 

! Alt p • c ‘'OPFN ( 'UNIT *3, NAME ■'■'ANCHOR, DAT '.TYPE •'NEW'" 

‘ A??P C ♦ , ACCESS* 'DIRECT '.FORM* 'UNFORMATTED', RECORDS! ZE-ANCREC 

[ 'AB3P 'CC"""* ,INIT7AI.SIZE*1RAP 

aBAP ' C. maxREC ■ R7P, ASSOCIATEVaRIABLE ■ INOX I 

r~ ~a BSP c 7”77’"ZZ 777 7 7 7 

' u jMTi c ' 

7' a B 7 p c ; Skip "To' s econo file " " 

A BSP '■■'C " ■ ... - - ... 

1 AB9P •---c--y<r "'CONTINUE'"' ■ ■■' 

" Aa;Vp ^CA'L'r s'K lP'"'(0NIT"2, IC)"'" _ '7 ' 

I" "AAip c "7' 7 7Z Z Z 

■ AA?p - CALL"'SK1P (UNIT,T,ICJ ■■ ' 

4ABP — ~'TF nC.NE.EnF) STOP ' !EXPECTE0'"EN0 OF' FILE ^ 

""AAAP ""'C 

f AANp ■■■ C'"' 

AAEP C4l'l''RFA0""(UNIT,LBF,532S,NAR, IC) ' 

[ a ATP C"' 

" AA8P IF {NAP.NE, 5328 1 STOP 'UNKNOWN FORMAT' 


! 
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SEO 

CZCSMBG.rTN ?6-PCB-Rl 12144107 PAGE 9 



It.... 



. 1 


C 




as(?0 

c 

CONSTANTS FOP HEADER FILE 


r as 1 0 

c 




as?0 


NE«198« 

1 NUMBER OF OUTPUT ELEMENTS 

1 osTa 


NL «970 

1 NUMBER OF OUTPUT LINES 

OSU0 


MHPo8 

1 NUMBER OF 8VTES/PIXEL 

[ 4SS0 


TPN«5 

1 PROJECTION NUMBER 

as(5{ii 


NPPI0 

1 NUMBER OP PROJECTION PARAMETERS 


[ «srd 


IFC«0 

1 element offset 


<198(1 


II O«0 

1 LINE OFFSET 

r as90 


NOWL«8 

1 number of overlays 


<lfi!»0 


0X»H?9, 

1 SAMPLING INTERVAL (METERS) ACROSS 

TRACK 

! a M d 


OT»805. 

I SAMPLING interval (METERS) ALONG 

TRACK 


c 



1 4 8 50 


NHaNE *NBP 



<1840 

c 




r a 6 S 0 


IF (LOAOl CALL IMOPtNn,ST,,lDR,NB,NL,'M») 


<1880 


OPFN ( UNIT « 4, NAHE • ♦ C SP. 0 A T * , T V PE • * NEW * ) 


1 <1870 

c 




<|8B0 


GO TO 99401 



1 <18 90 

c 




<I7(?0 

c * * * * • 




r <1710 

c 




O7?0 

c 




r a7T0 

c 

LOAD HEAOEP FfLE 



47(10 

c 




r 4 7S0 

99800 

CONTINUE 



4780 

c 




!■ 4770 


CALL GHFN fST.HFN) 



4780 


CAI.L ASSIGN (1 ,HFN) 


r 4 7 90 

■ c 




4800 

c 




r 4 810 


CALL M0ve{LRP(331.1,C00R0Sm. 1 . 0 , 0 ) 


4 8?0 

c 




[ 4«.50 

c 


• 


4 8 40 


SWAP « cnopnsc?) 



1 4SS0 


Ct.'ownsi?) a coawos(i) 


4 880 


coowosm ■ Swap 



I 4 8 70 

c 




4880 


SWAP B cnowos(4) 



i 4890 


COOROSfa) » C00P05(3) 


4900 


cnoRQSfS) B Swap 



1 4910 

c 




49a0 

c 




i 49 50 


LATLONCn * LATLONfl) • 9000 


4940 


IF (LATLONfPl.r.T, 18000) LATL0NI2) ■ - ( 38000 - LATL0NC2) ) 

\ 4980 


LT B FLOATILATLONfin /100 


4980 


LN » FLOAT a.ATLONf2) ) / 100 


( 4970 

. r. 




4980 

c 




! 4990 


wRlTFfi,i50) tear,oay;ht 



150 

FORMAT(lX,I4,':*,I3,?)i,30Al,7X) 

- 

i 9010 

r 




9000 

c 




i 90.50 


.WRITE (1,1801 NE.Nl 

,NHP,1PN,NPP,IE0, ILO.NOVL 


9040 

1 80 

FORMAT {819, 40X) 
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SEO CZCSMRC.PTN P6-PE0-8I 12l44»nr P*GE It) 


i 1... ...... I 1 .1 i 

C 

C 

'w»rfp n'.'i Ti, or 

sVhfi' (7W FORMAT_(«Pi2V>>V3PX) ” . •* 2 \u 1 

T' '5093 c 

■ 517.3 'C ' _ - ^ - - : 

[ 5113 ~wpiTr”nvi'70r tH,”i72,2^ ' " • * : : 

■ C“ "■ - " 

1 51T3 C '■ 

5 ia 0 wQlir ~ ' 

[ "SI 53 "TBr FORMAT' (1615) 7 T - ^ - - • - ^ 

■' '5(63 ""C '■ ' ^3" 

i ■■ stToi- "‘'call' CL bst""(t ) — — : " 

' '■" 51 8 »■■”"■ c‘ r'~"i • - _ r_ 

1 5193 .'GO "TO R96'0r " '■ * 

S?33 C""' ' 


5??3 C 


r 50 50 C" 

" 500 3 ■' C '■ {. 0 A 0 ' N A 01 H ri L E"" 

f-'S?5fl"”C 

" ■ 5 ? 6 3 99303 "C'OAj T T NU F 

p- 5373 C - 

50>13 "■ t ■ - 


r"'5093 ■'S'WAF"ir LBFTit) ^ ■ ■■ 

' 5533 TeFCirrvtPf'dB) “ ■ 

I " 5 T 1 '3 L'B F ( |'8 y-ii ■■ 5 * p - ^ : 

5i03 '"'~r ■" 

I ‘SWAP i"TPF(19r ' ■“ 

' 5 y «3 ■ “LarnR) '* 'lhf'cpo') ” ' ■' “ ■" ' 

r”5J50- L'B'F(i»'0)"A"'Sw'AP ■ ■ ■ 

'■' 5563 C PP ' " 

j 5)73 WRTTF (4 , !53)'''vrAR',0AVVHT "*■ 

5)53 ■ RTlMF’ * FtnAtC IMRBtf. V /■■"6'3303'.‘V6TBA3F~ ' ' r"URITS""A MXnuTES ’ 

(■“ 5393 R1.1ORK i" RTIMF '♦ "( FLOATC'SPAN ) /'"60333. ) J END OF' IMAGE " 

5033 "'PJ WPT IE (<l,905V RTIME, RWORK " f "TIME IMAGE STARTvS, ENDS 

r 5«]3 DO 93 B WHCG9P » V, « 

‘ 5A03 C CAlCULAIE TIME OF Ml NOTE ' MARK’ REL AT I VE TO GTBASE , UNITS '• MINUTES 

! 59)3 RTTMF » float s NAOTTMC WHCGRP )) - GTBASE 

- 5(1071 9?3~ WRITE (A, 905) NAOLAT (WHCGRP) .NADLON(WHCGRP) .NAOALT (RHCGRP) 

r'"5g53"'‘ * ,RT1ME 

" ' 5(163 905 ■■ FORMATf 'pF l?,A 'V’U? . F12.S ) - 

I 5073 — - WRITF (U .937)" TLTVAl. ~ 1"1' 

‘5AR3 907~ format (Fi a, 6) 

[ 5 A 93 Z~ : " ■ ■ " 

5533 c ' ■ 

r “ 5 5 1 3 G OT'TO -9 9 33 1 ^ ' 

5500 C" 

C 5533 ■■■""■ CAi'iTJA'iVi'ir^V* *.VA *iAi'***iVi***A¥A¥*** A** 

55 4 3 'C 

[ 5553 '”c ’ ■""< 'V""' 

“■ 5563 ~ C — ~ EV TRACT THREE BYTES iNO FLIP TO CONFORM TO DEC CONVENTIONS ” 

[-■ 557 ;, c" : “ 

'-- 5503 - 99733 CONTINUE ■■■ 


l'V '5593 C 

5633 COOROS CO)”* ""300 
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...... I..., 1 

SSIP COOBOSfSl • LPI f OFFSET ) 

Sb?(» COORDSf?) ■ LBF( OFFSET ♦ 1 _ 

I 5bJ? ■ COOBOSm '■ lBF (' OFFSET V 2 ) 

" 5fiO« c tFSt >09 ' neg'ati Ve'slGN ' ' 

|_'SbS!l' L.«f f .ANO, '"?B0 

S^b3 IF (H09K,NE,tJ ) COOHOS(a)_» "37T 

i ihie e ' • 

' ■5bSir T'61s97si'r ” ' _ _ _ J ~ _ I 

1 SbSfl"'" c 

STpnJ J C ___2 ^ --- ^ . 


5/214 C»**. .****.*.*•*««. **..**•* 

r 5»50 ‘1'C 'jN'r'LUOE *C2CSM9G,I^C 

5 / A 0' 'r iiiTiiVj, • • * . . t't * *"*' 


S7M1 c 

7 7 ?L H1__I iNO 


r 
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COASTAL ZONE COLOR SCANNER GROUND CONTROL PROGRAM 



SEQ 

C.PTN 

?6>PEB«6| 12I47M1 PAGE 1 








10 

C 




?0 

C 

COtSTAL ZONE COUOP SCANNER GROUND CONTROL » CZ&C 



40 

C 

. • • 

, 



C 

AnOPTEO ON OCTOBER 8, 1980 BV GLENN OAWIS PRON LCiCi' 



S0 

r 

REV 2,0 R OVE 30 OCT 80 TO INCORPORATE CONICAL SCAN 



60 

c 

REV 2,1 R DYE 11 DEC 80 TO USE LAGRANGE INTERPOLATION 


70 

c 




60 

c 





c 




1 00 

. .J,. _ 




1 10 

c 




1?0 

c 




1 ?0 

c 

**•* WAHNINGi TCHANGES IN MOUEl SHOULO ALSO BE HADE IN 

CZMP.FTN 


1 U0 

c 




1 S0 


REAL*0 LT,LN,LTLP,LNLP,LT0,LN0,LTSM,LNSH,S1D(4) 



1 60 


DIMENSION lYfl00), IV(100), PLT(10O), PLNfl00), EN(100), ££(100) 


1 70 


niMENSlON TALE(O), PALE(O,100), JOL(390), PALT(100) 



1 6 0 


■ DINENSION' A (10) 7 X(4) , X9(0), soy(4), R(0) ' ■ 



1 90 


dimension 6(15), Y(51, YB(5), S0V(5), P(5)" 



300' 


DIMENSION TT(0),TLT(0),TLN(«),TALT(0) 



?10 


■ INTEGER*? IPT (100), 1W( 100) ,1GN( 100), 



??0 


INTEGER*? AN,AE,4Y,BPE,0NE 

1 


P J0 


CORICaL*! TM(8),04{9) 



?a0 


DATA R£,RP,nBH/63T6200,, 6356800, ,57, 2957799/ 



i?S0 


data aN, ae, ay/ 'N* , 'E* , ' y*/ 



?60 


■'DATA 'SF/ 1000000;/ 



?70 


DAlA"HEV/2;i/ ■ ■ 



?60 


LF(IP,NIP, I, J)tJ-I*NlP*l-((iP-l*l)*(IP-I*2))/2’ 



?90 


ONE « 10 



T00 





T10 


white (5,20) REV 



1?0 

?ei ■ 

for mat * ( ♦0CO4STAir ZONE" COLOR SCANNER' GROUND CONTROL* ' 

REV » 


T10 


♦ ,F3,1,//) 



3«0 

30 

format (?F10,0) 



3S0 


RES»Rt*RE 



160 


EC»RE/HP 



170 


ECS*EC*EC 



360 


■ call* ASSlGKf (6, •DBICZCG,1ST») 



390 


LP«5 



000 

c 




ai0 


IALL ■ 1568 



O?0 


LLO • TALL / ? 



030 

c 




a «0 


t'AlL AS8iGr«i(0, *CSP,T5A7*) 



050 


READ (4,40) SIO 



060 

00 

F0RMATf4A0) ' 



070 


“ WRITE(5, 40)810 



060 


■ -READ (4,50)10 



O90 

50 

FnfiMAT(?Fl2,6,Fi2,0,F12,-6) '* ' * * 



500 


■ 00 60 1*1 ,4 ■■■ *“ 1 



Si 0 


* REAr)(4,50)TLT(l) ,TLfl(n,TALT(nVTT(n* 

• 


S?0 

"60 

WRITE (5,50) TLT(l),TLN(I),TALT(I),TT(n 



S10' 


BEAnf4,50)TILT'* * * 



500 


"CALL CLOSE (4) * * — ' 



550 


TSBNINTfTILT)/(2,'*0PR) ' 



560 


FF*45,/0PR '■ ' - - 
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SfcO 

C.FTN 

9b>FeB-ei laianii page s 


* • • • • 


STB 

C 


SB? 

c 


r 5«>p 


call assign (2, •OBIOSL.001*! 

(>(«(« 


OEFINt FILE 2(100, 18, U, INOlf) ' 

(lid 


CAIL ASSIGN (3,»0B|MCP.riAT») 



IG»1 

63^ 


mi 

bU? 


L«*l 

bSP) 

c 


bb^ 

c 


b7P 


N«0 

b80 

lb0 

CONTINUE 

bO0 


I«N*1 

700 


READ (3, l70,ENn«?00) M,PLT(l),PLN(I),ALT,(PALF(J,J1,J«l,ai 

7t7l 

170 

format n3,lX,2Fl?,b,F12.1,lK,4AA} 

7?0 


IF (M.fcO,01 GO TO lb0 

730 


I»NA1 

7a0 


RF.AO (2*M) IG,NO,TX(I),iy(l) 

7S0 


IF nG.EO.0) GO TO lb0 

7<50 


NsNbl 

7 70 


IPT(N1«M 

7M0 


IW(N)«1 

790 


1GN(N1«IG«1 

fl00 

C 


P10 

C 


«?0 


GO TO lb0 

m<n 

C 


««0 

?00 

CALL CLOSE (3) 

BS0 


WRITE (5,210) 

'Sb0 

?10 

F0RM4T (‘SWITCH AND ROLL DEGREE ? *) 

670 

' 

READ (5,?20) lOG 

860 

??0 

format (I5,?F10.0) 

890 


IF (IOG.LT.I.OR.IOG.GT.3) I0G»1 

900 

c 


910 

c 

• 

9?0 

?30 

CONTINUE 

930 


KHITE(5,a5n 

9a0 

?51 

FORMAT ( 'SI.NSH,LTSM,VM, tilt 7 *) 

9S0 


RtAO (5,?53HNSH,LT5H,VM,T8 

9b0 

?5? 

FopM4t(<IF10,0) 

970 


TS»TS/OPR 

980 


00 9(10 JI<1 , 10 

990 


KbI 

1 000 


no 2(10 III ,4 

1010 


yem« 0 . 

10?0 


IF (ji.Eo.n R(n« 0 . 

1030 


DO ?(I0 J«I,4 

1000 


A(K)*0, 

1 050 

PO0 

K«K*1 

1080 


KbI 

1070 


bd 250 ibi,5 r 

1080 


V8(h'«0, 

1090 


IF (JI.EQ.l) P(n»0. 

1 1 00 


on 2S0 J*l,5 

1110 


P (K 1 B0. 

1 1?0 

?50 

KbK*1 
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SEO 

C.PTN 


?b>PEB*81 

12I4TI11 

PACE 3 




I 

1 • • 1 • • • • 1 


ii 1.. 





1130 SSN«PI, 

II 0 3 S 5 1 ■0 , 


! 


sw» 0 . 




1 


CTEWB0, 


, 


r 1170 


CTEWP«0. 




1 1 pp> 


CTNS»0. 




! 1 i9«i 

c 

' 




1 


no P8?l 1CP«|,N 




1?10 


St 1 FLOAT ( iV(XCP) ji - 485. 1 N0RMALI7E 

ALONG 

TRACK 

dimension 

l??7i 


ST « ST /485, 1 normalize 

ALONG 

TRACK 

DIMENSION 

1?50 

c 





1 ?<I0 


t ■ ( 9 70 . - 1 V (I CP ) ) * i ?8 , / (970 , *60. ) *70 




)?Sil 


CALL LR»(TT.TLt,t,LTO) 




1 


CALL LGWrTTjTLN.T.LNl)) 




|P70 


CALL LGKrtT,TALT,T,ALTJ 




l?«0 

c 

WBITe(5,50ILTO,LNI), ALT.FLOAtdVilCP) j 




1 ?9CI 

c 





1 JPI0 


LT * LTO / DPR 




nn 


In ■ LNO / DPR ■“ 




|3??( 


LN»LN*LN3H 





i 1 3S0 ' ITtUT + LtSH 


— 

1 Tu0 


■SLTsSIN(Ln 


1 ?S0 “ 


CLT«CaS(LT) 

; 


1 3M 


SM«,15839/CLT 



1 370 


CH*.S0PT(1.-SH«SH) 



i 300 


HD*ATAN?(SH,CH) 



'1 390““ 


Hn«HO*VW 



1400 


Hno«HO*i)PR 



1410 

C 

* 



1 a?0 


4L«AT4N(SLf /(EC*CLTt j 



lil30 ^ 


~ R*L*RE*COS(ALj+ALT*Ci;t 



1440 


X0«HAL*COSfLN) 



1 4*50 


V0*H4U«3IN(LN1 



1 460 


Z0*RP«SIN(AU*ALT*SLT 



1 4 70 


C»PAL»Ra|, ♦70*Z0*ecS-RES 



1 4 80 


GCLT»atan(Z0/hau 



1490 . . 

C 




1S00 

c 




1510 


PL» ,00 



15P0 


RL *PL ♦ ( ( f R ( 4 ) • 8 1 AP ( 3 n • STtR (2) ) * ST*R { 1 ) ) ? ALT 



1S30 


~ PT* ,0 



1S40 


PTiPT*(((Pt5t*ST*P(4))*ST*P(3))*ST+P(n)/ALT 



1550' ■ 


ho»md-P(3) 



1560 

C 




1 S70 


pNiiKnCP)-LLT) 



1 580 


• Tri»PN*VM/ 1000000, 



■'1590 


*1«0, 



1 600 


T1»0,' 



1610 


71*aLT 



1 6?0 


call SARtm ,71 ,T3,K?,Z2) 



16 30 


>?*YJ * 



1640 


CALL SARrv3, K?,-TM, 1^3, K3) 

' 

. 

1650 ■ 


73SZ2 



1 660 


CALL S4R(X3,Z3,FF,X4,Z4) : ' “ ” 



1670 


V4*Y3 



1650 


"■ y4«-Y4 IMIRROR REFLECTION 
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SEO 

C.FTN 

?6«PEB«ai 13t47tll PAGE 

a 


t • • • 1 



itsqo 


CALL SiR(xa,2a,.FP,V5>ZS) 


1 TCTff 


¥5«va 


[ ITllpi 


CALL SARrV5,*5.TH,V6,X61 


1 


78125 


i Tiff) 


CALL SAKM(8.?8..TS,)f7,Zn 


1 Tapi 


V7«V8 


1 7SB 


CALL SARfX7, V7,-»Lf *3|731 


1 76PI 


Z3»Z7 


177(* 


call sab (73.73.PT.Z3.X3T 


1 7 8« 


CALL SAB (TS.ZS.HO.Va.Za) 


t 790 


CALL SAB (Z9, X3.GCLT.Z5.X5T 


1 800 


CALL SAB (Va.XS.LN.Yfi.Xfel 


1810 

C 



1 8?0 


AQaX8*X84V8*V84Z5*ZS*ECS 


1 8A0 


80*3. A (X6AX0«Y6*y0AZ5*Z0*ECS)/AQ 


I 890 


Ta.5* (•BO-SQBTlBOaBQ-a.aC/AO)) 


1 8S0 


XEaX6*rAX0 


1 880 


YFaY8aTAY0 


1870 


ZEaZ5*T+70 


1 880 


L IL^aA r AN(ZFafcCS/SOBT(XEaxe*YE*YE)l 


1890 

* 

LNLPaOP9*ATAN3(YF.XE) 


1 900 


LTl PalJBBaLTLP 


1910 


ELT*LTLP-PLT( ICP) 

j 

1920 


ELNaLNLP-PLNdCPI 


1 9.80 


WTalwdCPl 


1990 


FN!aeLT«9F./l1PH 


1 950 


EE laeLNABeaCLT/OPR 


1 980 


CAI.L SAB (ENI.etl.HO.PTE.RLE) 


• 1 970 


RLEaBLEaCOS(TM) 


1980 


XfDaHLE 


1 990 


X(3)aST 


?000 


X(3laST»8T 


P0] 0 


X(9)aSTa8lAST 




Kal 


?050 


00 380 lal .a 

V . 

P090 


X8d)*XBdl*Xd)aWT 


?05R 


00 380 Jal.a 


?080 


A(K)aA(K)AXd1aX(J)aUT 


?070 

• 380 

Ka« + l 


?080 

C 



?090 


ymepTE 


? 1 00 


Y f 3) aST 


?1 10 


¥{91aSTaST 


?1?0 


Y(S)ssTaSTaST 


?1 30 


Y(3)aTHaALT 


2190 

c 



3150 


Kal 


3160 


DO 370 1*1,5 


. 3170 


YBmaVBmAVdJaWT 


3180 


00 370 Jal,5 


31 90 


P{K)8B(Kl+YdJ*V(J)*WH 


3300 

370 

Kaitl ~ 

- 

3310 


RwaSWAliiT 


3330 


SSNsSSNAENTalNIaWT 


3330 


SSFaSSEAEEIaEEIaWT 


3390 


EM(ICP)aENI 
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SEO 

C.PTN 

36-FEB-Ht 13I47IH PAGE S 


if «... 


??S0 


Eencp)»efci 



CTEW«CTEW»TH.RLE 

\ ??70 


CTEW3*CTEW3^TM*TM*RLE 



CTNSiCTNS^TM»TM*PTE . ■ • * " 

L . 

380 

CON fl NIUE . s^--* 


C 


1 P510 

C 

‘ 

?S?9t 


CTEWaCTEW/SW 

i 


C rEM?iCTEW3/8U 

?sa 0 


CTNS'CTNS/Srt 

r PTSfi 


Kil 

?«60 


no 390 !«!,« 

r ?T70 


DO 390 Jal ,4 

?T »>0 


AfK ) sA (K)-KP (tl *>(»f Ji/SW 


390 ■■ 

KiK + 1 

?a 00 


Ktl 

f ?«lfi 


00 300 lal,a 

?«?0 


)(H(n»><0rn/8w 

r 


00 300 Jat.U 

?aa 0 ' 


AfK5iA(K)/(SW-l,) 

•■" 'pasi* 

300 

KaK + 1 



no'3i0TaiV4 * "■ r" 

"" ?a70 

■ 3 1 e 

SOxrnaSOWTfAflPfiT, 10,171)1) “■ ■ ■ r 



Ka 1 

P't'jq 


no 330 la 1,4 , 

?S 00 


on 330 Jal,4 

■" PST 0 


A(K)iA7Kl/(30r(ITa$n)<(J))' ■■ ■ ' " ■“ ■ 

?s ?0 

*330' 

KaK.) 

‘ ?SS0 


00 330 )a3,itI0G 

3 SU 0 


IK i:*(L>(a,i0,i,n),LT,j,e-5) so to T40 

?S50 

" 530' '■ 

CALL step (A, 10,4,1) 

?560 

"c 


3S70 


vniaxem 

?S80 


00 340 Ia3,l^I0G 

?S «0 


*maSnxfn**(LK(4,10,!,l))/SDX(n ' 


340 

xf nsx(n-xm*xH(i) 

?61 0 


DO JS 0 laljl.IOG 

?6?0 

3S0 

»(I)aRn) 4 XCI) 

?fe ?0 

c 


?f.un 


Ka 1 

36*50 


no 360 I a 1,5 

?660 


on J (,0 Jal,5 

3670 


flfk).B(K)-yB(n*Y0(J)/3W 

3hfl0 

360 

KaKtl 

3690 


Kat 

3700 


On '370 Tar;5 - ■■■ ■■ 

37 1 0 


VBmaV0(l)/9W ■■ ■ ■ 

37?0 


no 370 J a 1,5 

■■■-P7 30' 


P(K)s0(H)/fSW-J.) 

' 37U0 

■ 370 

KaK.) 

'■'3750 


00 330 Tai,5" 

' 3760 

'-380' ■ 

soy(nasoRT(0aK(5,i5,'i,i))) 

' 3770 


K ® J 

■ 3790 


no 390 I a 1,5 

■ 3790 


DO 390 Jal ,5 ■■■■ " •. ■ 

3800 


HfKVaBfK)/(SnY(nASDV(J)) 



“ 
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PfllB 

390 

K«K41 


?8?« 


on 000 I«2,2«10G 


?8T!? 


IF (8aF(5,i5,I,lh,LT,l,e-5) go to TO0 


?9II0 

000 

CALL STEP fB, 15,5, I) 



C 





Y n 1 ■VBM 1 


?»70 


00 0)0 la0,2«lOG 


?8H0 


vm.SOYf n»fl(LF (5.15.1.1) l/SDYfll 


?890 

0] 0 

Yfi )«vm.Ynj*vB(i) 


?<»5in 


00 0^0 Is),04TDG 


?9t0 

fl?0 

piti)»p(n*v(n 


?<»?0 


S5N»SSN/SW 


p'lia 


SSEaSSfc/SW 


?9a0 


SON»SOHT fSSM) : 




SOtaSO«T(SSF.1 


pQbB 


W9ITF. (5,O30)SON,8OE 


?970 

030 

FrjPM*T(0F15.1) 


?9H0 

000 

CONTINUE 


?990 

C 


r ■ 

5000 

3010‘ 

■ — 

WHITC(5,001 )CTEW,CTEW2,CTNS 
FOBMAT(3ei5.5) 


30?B 
50 5(i‘ 

C 

005“' 

CALL DATE (6A) 


50<I0 


CALL time (TM) 


5050 


W»irE (LP,O501 0EV,O4,TM 


3060 

050 

format (M coastal tone ground control points REV '.F3.1.5V.9A1.2X 


3070 


i,8Ai./7i 


3080 


WRITF (LP.O60) SIO 


5090 

060 

format f* SCENE 10 *,0A0,/) 


31 00 


WRITE (LP, 0651 LNSH,VM,TS»0PR 


3110 

• 065 

F0BMAT(3F1?.5) 


3I?0 


WRITE (LP,O701 SOE,SON 


3130 

070 

format {' RMS ERRORS EAST ».F7.I.* NORTH •,F7,1,*’ (METERS)») 


3140 


WRITE (LP,a601 P(n,P(3),P(fl).P(5) 


31^0 

080 

format f'0PTTCHl',OF10,ll 


3160 


WRItE (LP.O901 R 


3170 

090 

format (' ROLL I'.OFIO.!) 


31 80 


WRITE aP,500) P(2)»1.E6 


3190 

500 

format (' YAW l'F10.n 


3?00 


WRITE (LP.5101 


3?10 

510 

format (//' POINT FILE WEIGHT EAST NORTH',/) 


3P?0 


no 530 1 8 1 , N 


3?30 


NSOE«nINT(EE(I)/SOE) 


3P40 


NS0N»NINT (ENd) /SON) 

1 

3?50 


WRITE (LP,5?0) IPT(n,IGN{n,IW(n,EE(I),NSOE,EN(I),NSON, (PALFfJ, I 


3?60 


1),J»1,0) 


3?70 

5?0 

format (IO,I5,I6,FIO,1,IO,F12,1,T«,OX,OA9) 


3?80 

530 

CONTINUE 


3?90 


IF (LP.E'0.6)G0 TO 610 


3300 


IPWF«0 


3310 


write: (5,500) ! 


3 3?0 

500 

format ('SALTER P0INT\WEIGHT 7 ') 


3 3 30 

550 

REAO (5.560) IP,ITW,FFF 


3300 

560 

format (?I10,F10,0) 


3350 


IF (IP. 60,0) GO TO 590 


3360 


tTWiTAHS(ITW) 


70 




ORIGINAL PAGE 19 
OF POOR QUALITY 

APPLICATIONS DIVISION 


seo 

C.FTN 

06-FEB-ai 12I47M1 PAGE 7 




3370 


IF (tTW.GT.lB) 1TW«10 

3380 


00 5T0 I»1,N 

r 


IF (TP.fQ.IPTdn TH(n«ITW 

3a00 

570 

CONTINUE . . • • ■ 

S«M0 


iPw8 •! ■ . '• ' 

3«?0 


WRITE (5,9801 

3030 

SS'B 

FOR'MAt I'SANdtHER 7 ^ 

3040 


GO TO 990 

3080 

990 

IF (IPWF,NF,0) GO TO 330 

3060 


WRITE (5,6001 

3070 

600 

FORMAT (•# OUTPUT ON LINE PRINTER t ♦) 

3080 


RE4U (5,6091 LPA 

3090 

605 

F0RM4T(A11 

3S00 


IE (LPA.NE.AvJ go id 610 

3910 


LP«6 

39P0 


GO TO 005 

3830 

6t0 

WRITE (5,6?01 

3900 

6?0 

format^’ I' sPRiNrER Plots t m 

5990 


'REAiy"(5,6051 LPA'. 

■'3960 


TF (uPAVNEiAyi GO"TO'720 ■ ' 

— 3510 

C 


■ 3980 


■WRITE (5, 630) ' “ ' ■' ” ■' " 

— 3990 

^6'30 

"Format ('Smeters per oivtsion't »)' " " ' i 

3600 


WEAO (5,301 PPS 

■■■'3610 


IF (PPS,F(J,0,i PP8«P00, ■ ■ ■ 

■-■ 36 ?0 


PP8FtPP5/l0, 

3630 


WRITE ■ (6;6a0r'PP3'"" ■■ — 

3600 

'6‘ii‘0 

• format ('tPlXEl. v;8 " errors *',T 28',*E AST*, ■T88r»NORTH*7Fl0,0r*METERS7" 

'■ '3690 


lOIV'/l 

3660 


00 670 1*1,60 

— 36 7 0 

' 

'LL*CT-n*33 /' " ’ '■■■■“ 

3680 


LU»I*33 

3690 


on 660 J*f, M 

3700 


IF riw(jl,eo.01 GO to 660 

r 3710 


IF (IX(J1 .LE.LLI go to 660 

37?0 


IF (lx(J),GT,LU) GO TO 660 

r 37 30 


leE«Ft(Jl/PP9F+30, 

37 00 


IEN*EN(J1/PPSE>90, 

["■ 3790 


"IF-"(lEe,Lt,Pl IEE*3 " ■ 

" 3760 


IF (IEE.GT.601 166*60 

; "3770 


IF (I6N,l.T.6n IEN*61 

3780 


IE (ItN.GT.tPBl IEN*130 

3790 


"WRITE (6,6501 1PT(J),IPT(J) 

1800 

650 

FORMAT (•♦',T«IEE>,1?,T<IEN»,I2) 

! 3810 

660~ 

CONTINUE 


6TB' WRITE' rb','6«B) ' ■“ 

r~"? " 30 6 P 0~~P 0 R M 4 T~ ( 5 ( R X 1 * T , 1 0 X , 5 { 9 X I • ) ) “ 

l^HlTE't 6,690) " 

r“3950'~~69'0" FORMAT" (MIT ME"'V,S,- ERRORS •VT28;'»E»ST*7T88pN0RTH*;') 

— '3860 on 71B' Iat,feB — 

r 3870 LL « 1 1-1 ) * 1 T """■ ' 


3880 

LiJ< 

>LL*17 





- 

' 3890 

no 

700 J* 

1,N 





■■ ■ 3900 

'■ IF 

(IW(J1 

,eo 

.0) 

GO 

TO 

700 

■ 3910 

IF 

nv(Ji 

.LE 

.LU 

GO 

TO 

700 " - - 

3930 

IF 

(IKJl 

.TiT 

.LUl 

GO 

TO 

700 "■ ■ " "■ ■ 
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i.... I 1.. i t 

iee«et(j)/ppsr*30 

_ieM»pNf J1/PPSF^9P, 

J?_riEe,Ll,?) IEE«3 

IP (IEt;,i?T,»i0ii l€b«60 _ _ _ „ - _L 

"J_I_I1rP nEN.LT,6l) lENifcl — - "" - 

IP tl£N,r,T,iPW) IEN« 1 ?« __ ^ 

___W«ITF (fc,hSP)_IPTIJ),IPtJiL) : 


flMCTP T 00 CONTINUE __ _ 

■ acii(?i T'ta'^' wPiTE (tijbaoj I IILIZI—ZI™ 

~ I ^ , 7 p ( I J 

■ ae>l3P" ' ^WPITP (<l,730) R ' 

■floom ■ WRITE (fl.TJBI P(?)ii,EB ' _ 

«PS0 '^3«_l_fOR'^*T (BMP, 31 

aPbP ^ _ ■ ’close (IJNI T»fe, 0!SP08EV»PPlNT*) _I ,Z_ IZ . ' Z_..'' 

" «P70 ZSTOP _ Zl_ Z'_ Z_ Z_ 

aosp faa'- wwtte. f'5,7S0)’'“ 7 ” _____ 

^aPRP . - ■ ^ 0 RM|,T ' i ' 0CP„P>f i„NOT BUITABIlE j'MJEM ~ 

_ aiP 0 ' ~ _ ZsTOP' " _ _ _ 

' at 10 " '• ' ENO T' Z Z ' 1 

“ aipp “■ SUPROUf INE 'lCR( ’ v f PV , i'j , 'iiRESLT) _ _ _ 

■^ats'p’Z'C - - -- “■■■■■_'_ 

a I ^ei" t ' iizzi zi" ' z _ I - zziz” z_z zi r i_ zz 

■“aifep' _c pn'vau _ ' ” ^Z- Z _ 

~ at70~" c XZZ" _ Z I~ZZZ iNVAtsL- .ZZ ^ 

^aip 0 c J _ __*Rf 8 ».f • _l _ FNVA|.a - _ __ 

_a 1 90'Z’c ' iiS ^Zi I lFNy*Xi.?. I i-.,.. 

a ?00 c j I i, I ; 


a ??0 

c 

« «» « * * 

a?S0 

' c 

VALl - XINTRP • VALP <• VAL3 - VALA • 

aaap 

c 


a?s0 

c 

GIVEN POUR VALUES ( VALl -> a OP ARRAY V ) AND 

uab0 

c 

GIVEN POUR corresponding VALUES ( PNVAUt -» a OR ARRAY PV ) THAT 

a?70 

c 

are a punction op the pohmer values. 

a ?80 

c 

ANO GIVEN A value VINTRP ( XI ) THAT IS ADJACENT TO VALl -> « 

a?90 

c 

THEN THIS routine RETURNS A VALUE, XRESLT, THAT IS THE LAGRANGE 

aS00 

c 

interpolation op xintrp about fnvah -> a. 

aj»0 

c 


a 3?0 

c 


as30 

c 


a3a0 


real via) , Fvia) , xreslt , xi, PART(a) 

aS50 


REAL A.B.C.D 

a360 

c 


a370 


A"”. XI ■ vin 

a380 


R • XI - VI?) 

r a 390 


C P XI • V(3) 

aa00 


0 p XI - Via) 

aa 1 0 

0 


aa ?0 

0 

type S00,A,e,c,o 

aa 30 

0 

500 PORHATiaPia.B 1 

aaa 0 

0 

type 500, V ' 

aas0 

D 

TYPE 500, PV 

OObPI 

0 

TYPE 500, XI 

aa70 

c 


aap0 


PARI in P I B A C p 0 1 

- 
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i... 



. 1 ,. 1 .. 


asc^0 


♦ / (( v(n • v(sn * Mh) - 

♦ * PV (1) 

V(31) * ( V(l) • 

v(on ) 


«S10 

C 





as?(>) 


P*RT(?) ■ { A * C * 0 ) 




os 


♦ / f ( v(2) - VM)) • ( V(2) - 

V(3)l. i .. 

v(«n ) 


OSltPl 


♦ * FV (21 



05S3 

C 

• 




OSbti 


PART(5) ■ ( A * H • 0 ) 




0ST(>l 


♦ / (( V(31 • Vnn * ( V(3) > 

V(2)) * t V(3) - 

v(on ) 


OS8i» 


♦ * FV (3) 




0S9H 

C 





06?tP 


PARt(0i ■ ( A * R * C ) 




Obl« 


♦ / (( vffl) - v(in « ( v(o) • 

V(2)5 * LV(fl) - 

V(3}) ) 


ubpn 


♦ * FV (01 




OhSPt 


VRfSLT • PART(11 ♦ PART(2) ♦ PART(3i 

♦ PART {9) 



0b0l>l 

C 





065B 

0 

TYPE ■>00, PART 




Obbffi 

D 

TYPE S00,XRf8LT 




ObTS 

' c - 





0880 

C 





0890 


■"■■■RETURN' : 




0700 


END 




0710 


'SimROiJTlNE'STEP (A|NtF,IP,KAYr" ■" " 


; 


O7?0 


OlMENSldN A(jSl 




0 7! SI' 


tF(U, J*T«Jil-IXYNiP+l*((lP«I*9l)AClP 




0700 


LkK*LF (KAY.KAYI 




0750 


H«0 




'0780 


00 70 I»1 , IP ' ■ 




07 70 


Oti 70 Jo I, IP 

' 



0780 


MCM«1 




0790 


- 1 K ~( I - (CA Y 1 " 10, 7 «V99~ ■ ■ ■■ 




0800 

10 

LTKsLF(I,kAy1 




0810 


GO ro 30 




OH?0 

?0 

LIk«LF(yay,I1 




0880 

S0 

IF (J*KAyi 00,70,50 




U8U0 

00 

LKJ'LF ( J,K*Y1 




O8S0 


GO TO 80 




0880 

50 

t.KJ«LF (FAY, J1 




0 8 70 

60 

A (h) *4 (Mi- A (L IK)* A (LKJ)/A(LKK) 




0880 

70 

CONTINUE 




0890 


DO 110 1*1, IP 




0 900 


IF (T,K4Y1 80,110,90 




0 910 

80 

LIK»LF (I.KAY) 




O9?0 


GO TO 100 




O'950 

90 

ClK»CF(K4Y,n 




O9U0 

100 

A(LIKi«A(LIK)/A(LKK) 




O9S0- 

li0‘ 

CONTiNlIf 




0980 


*(LKKl»-l ,/A(LKK) 




0970 


RETURN 




0 980 


END 




'r 


\ 
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C75MP.FTN ab-6EH.Sl 12I00I8B PAGE 1 

C 

COASTAL 70NE SCANNER MAPPING POLYNOMIALS 

c 

AOApfEO FROM MMP 9¥ GLENN OAVIS ON JANUARV 13. 1981 


c 


j <10 


logical AOVFL 

50 


Om.iHLE PRECISION Y.XN.A 



DOLirtLE PRECISION 5ME AN, STOEV , POSO 

70 

L 

90 

— 

DIMENSION V(OOI, Afsao), CtOO), FfOO), FE(Ot)) 

DIMENSION KMEANCOO), STDEV(OO), B(O0), 0(00), TOLEVtOOl, P(00), NI 
Tfnm0i, iNENfiicn 

[“ 100 


DIMENSION TT (0) ,TLT(0) ,TLN(0) ,TAlTt0) 

1 10 


integer stn(5 ) ,ay,rpf 

[ 1?0 


LOGICAL*! I T (101 ,MFN(181 ,TM(8) ,0A (9) 

1 50 

150 

[ IfcO 

1 7 0 

[ 180 
190 


real LT,LN,LTLP,LnLP,LT8C,LNSC,PT,RL,LTSM,LNSH,SF 


real LTP,LN0 

— 

DIMENSION PS(?) ,PnCH(5) , ROLL ( 0) , delta (2) 

common X, TN, A,C.F,Fe,L,NTGC,IP,NIP,FINC,FbUT,KAir,FLAG,KOEP,6F,tOL “ 

DATA RE, RP, DPR, SF/63T6200,, 6356800, ,57,2957795, 1000000,/ 

DATA HMVE/OHRMVE/.ENTR/OHENTR/ 


DATA AY/*y*/ 

1 ?00 

LMiK, j)o»jx-is*Nip*f-((ip-ix*n*np-ii(62n>2 

?10 


4DV6L«,TRUE. 

r ??o 


KBaS 

?50 


LPa3 

r ?ao 


CALL ASSIGN(3, »CZSMP.L8T») 

?90 


WRITE (kb,?0) 

r ?6o 

?70 



FORMAT I'OCOASTAL ZONE SCANNER MAPPING POLYNOMIALS REV 1.0'/ 
1 'SPROJECTION NUMBER ? ») 

|- . ?«o 


READ (KP.301 IPN.LR 

, ?on) 
r '"500 

30 

FORMAT (?nO) 


RESaRE*«E 

510 
r“ '5?o 


ECaWE/RP 


ECSaECaEC 

5 50 

C 


! 5u0 


lALL a 1968 

550 


LLO a IALL / 2 

r 5b0 

C 


570 


WRITE (8P.57) 

[ 580' 

37 

format i'simage title T ') 

390 


PEAH fKH.sn IT 

f 000 

31 

format MOAl) 

010 


CALL GMFN (IT, HEN) 

r 


CAIL ASSIGN (2,HFN) 

U 50 


RFAU (?,5?) sin 

1 0 00 

5? 

FORMAT (1Y.5A?) 

050 


read (2,33) NIE,NiL,N(5P,IIP,NIPP,IE0,lL0 

i 060 

33 

FORMAT (715) 

0 70 


READ (2,3<n nELT4(n,0ELTA(2),LTSC,LNSC 

[ 080 

50 

FORMAT (aE12.61 

090 


CALL CLOSE (2) 

1 500 


CALL GT(LTSC,LN8C, IPN;U0,V0) 

510 

C 

\ 

[ 5?0 


CALL A5SIGN(a, 'CSP.OAT') 

550 


READ (U.aOl SIO 

1 500 

00 

F0RMAT(5A?1 

550 


wRiie (5, aoisiD 

( 560 


REAIL (9.50)10 


f 
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SF.O 
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,FTN Pb-FtB.Rl 12144108 P4GE i 



* 1 • • • 



■57^ 

50 

F0«»«4T r2Fl?,6,M?.0,ri2,6) 

z: 

SP0 


pn 60 r*i , 4 

SPPI 


BEAb(4,50)TLtf n, TLNf I) ,t*tTm,TTfI5 

hPip 

60 

WRl7fc(5,50)TLTm,TLNm,TALT(n,TT(il . * ‘ 

z: 

hi 0 


Re4Of4,50)TIl.T 


6?0 


RFAD (4,7305 PITCHfl) ,0ITCH(3),PTTCH(45,PITCH(5) 


hS'’' 


RE»0 (4,7301 MOLL 


8ut1 


RL*0 ( 4,7301 PITr.H(2) 

j. 


750 

FOMM4Tf6F10.H 




CALL CLOSt (41 

r 

b7?l 


TS«NINT(TlLT)/f2,.0PRj . 


h«Ji 


FF «45. /PPM 

r ■ 

690 

■ C ' 



7 0 0 

c 

' 

! 

710 

7?0 



IP ■21 

"MIPi(iP*(IP-f|})>2 ■ 

r 

7 J0 


CALL ALtPPI(,45) 


7a0 

■■"C 


z: 

■750 


'CALL' assign (1, 'COEP.GeO'l " — ■ 


"'760 

■770 

- 105)-' 

“WRITf "(1 , 1001 SI0,IPN ■ ~ 

•FORMAT ' (SAP, 1101 


■ 7S0 

■ C 


[■■■■ 

790 


"CAu: OATH (D*y - - - 


900 


CALL tlMK (TM) 

1- 

810 

800 

"Ti 0 “ 

WRITE (1,1101 0A,7M 

■F0RMAT“(rx,9Ar,2)(,8Ari ' ~ 

[ 

"■*80 


■•■fc<RITE"'{LP, 1201 0A,TM ■■ * — - ■ ■ 



■'■C'"' - 

, 

z. 

890 
: 860 

1?0 _ 

format (»lCbA3TAL TONE SCAHNER MAPPING POlVnOMIaLS '9 A1 , 2X , 8 A1 /> ) 
"'WRITE '(LPVlSOl SID, IPN " 

f" 

870 

■■■■ 130 

"Format'’ ! ' SCEME'"IO ',5A2,'» ' "' PROJECTION NUMBER' , 'T«r’ 


880 

■■""r” ■ 


c 

890 


WR!TE(1("»,25H 


900 

251 

■FOhMAT('SLN5h,LT9M,VM,riLT’ T " " ’ 

Z' 

910 


REAO(I<0,25?1 LN8H,LTSH,VM,TS " ' 


9?0 

~ ?52 

FORMAT(4F10,0) 

L 

9?0 


TS « TS / OPR 


9a0 

■ C" 


j 

99:1 

■ c 



960 

C" 

IMIUALITE ACCUMULATORS AND MATRIX A 

r" 

L 

970 




980 


00 660 1TT»1 ,? 

l: 

■ 990 


M»0 


1 0j»0 


00 140 I»1,IP 


1010 


" XME AWnT ib,' ■ 


10?0 


“'00 140 J»l,IP 

i 

1030 


■ M s M ♦ 1 - - . . - - - 


1 0U0 

■ ra 0 ~ 

ArMi«0, "■ ■ ■■ 

z 

1053 

- f, 


1 060 


N»0 

z 

1070 

c 

• 


1080 


"00 "1 80 TNS 1 r~, 27 T' . - . - - - 

r 

L_ 

1090 


" TV "'» ■ (INS'-I 41" ' *■ 35^ ♦'"«85 ' " 

1 1 00 


ST ■ FLOAT ( IV 1 « 485, ' NORMALIZE ALONG TRACK DIMENSION 

i: 

1110 


"ST « ST /485,' 1' NORMALIZE ALONG TRACK DIMENSION 


1 1?0 

C ■ ' 


z; 



'• 


czz 

I 
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3EO C2SmP,FTN IPl4aiP6 PACE 5 

1 I.. ..I ..I ..I J J 

113a T»C»Ta,-!y)*l?H,/(970,*6«),)*T0 

ll«?l _ CALL LGRfTT,TI.T,T,LTO) 

[_il5a CAl.L_Li;P(TT,tLN,T,LNO) ■ 

llfeia _ CALL LGPrTT,TALT,T,ALty ■ ' “7“^ — 

■■ - ” - -jp atj,LT,3i uorn M? " 

T 'wRlTt(5,50»LTO,LNO,ALT,>L6Af nV) - - 

_ 1J2_ CONTI^NUt _ ' 


L. l?l'i _ J "LT ■ LTb’ / DPR 

__ i?pa -l’n * L*^b 7 OPR ■■ ■ ■ " 

[■ Tpjy' ln«l'n*lnsm ■ ■ 

Lt»LT»LTSM y ■ , 

1 iP'ia sltssin(ltii “" 7 I ’ 

'i?ha ’J"ci.T»cos(Ln '7“ T _ 1' I. *• ” ' 

L..J?7fl7 'll SHi,J5939/CLT 

'l?fla”‘ '"CH»-soRrn,-yH*SHV ■■■' ■ ■ - - 

L.ll?'’^ 1 ' HO«ATAN?fSH,CM) T 

1 Tpa ~ Hn»Ho*Yw "" 7 ” 

r 1510 II 1 _I 

lAppi' c7'7l 3 7 7 ■ 

i 133^ '■ ■■ AL«ATAA<f StT/7EC*CLf l ' ' 

--rsao 7 ' ■ RAL*Pfc*r, 0 S(AL)+Ai.T*cLt 

[' 1 55a " x 0 »RAL*cos rLNf ‘ 

’"15^P'7_ _ _ "¥0tRAL*SlN iLnY ' ' ' ’ ^ 

r''l570' ~ ■"■■■Z0«HP»3IN(AL')VArTA8Lt - 

1 ifiO C0 *ha”l»BAL47(»*20»EC8*RCS 

f7’ t5R0’_ GCL r«AtAN(70/RAL) ■ 

"i «00 7'c" 77 '"‘' 77 ^" • 

r'Voi^ 'c' ' 

_ia?pi’ _p[»' , 0 "ji _ __ I _Z_ IT _ 

[_'] a 5a7 IlPLiHI,.V'( TLRJILC i’p'syt^ (3f)isT*R0LL (2) ) PST'tPOLL (nj /ALT77_" 

1 U «0 PT« ,0 _ _ " ' 

r '~ia50 ' PT»Pr7(f'fprTC'H'(5Y*'s'TTpitCMf9)?«ST»P!TCH73 

lafep ITCH 7^1 / SF _ 7" _ 1" 

! 1 a 70711c Z! 11 . 1 1 1 'lz ir ' • ■ 

t ana' " " no i sai i'p'w » ~ 1 _ l’_ 1 ’ 11 

r ~ Y<i90 • iih:v~T»y' * 70 ♦'■pa'a * 

'“isa'a" ' PN - rUb” ” " 

r " l 5 1 0 _1_ 3 H * P N * V M / S F ' ' - — - 

7 VYi'ti ," ' ' ” ' ““ 

f~"l5 50”' yi»a," ■ - - ^ - 

'i5a0 ' “ 7i»ALT ’ ' - * 

L'T550 '■ 'c7CL''SA'R(Vi;z'f,rs', 7 FiY2^ ~ 

I5b0 T? = Tl" ' ■ ' ■ ■ 

r“l570 ' CALL SAPyv?,XP,«fH,V3;X3) ~ 

15»J0 Z5»Z? ' ’ 

7 1590 “ CArL SARfX5,25 7 FFV*a 7 z<iy - - - - 

If, 910 ¥««T5 ' — - _ 

[1 1 ^ ’ <’1_1 .: * a e - X iT ~7‘ 1 ‘fi'T RROR REFLECti ON ' 

i 6?0 r«LL SAR(X0,Z4,*FF,X5,Z5) - - -- - 

f “'1830 ~ y 8 »Ya~ 7 --•••• - - 

ibab _ __ y HC~sARrY'5,757fM7Yh,ia) _ ^ ^ 7 ” __ 7” ~ . 

1 1850 ' ~ ' Zb»Z5 77 ^ ■ 7 — 

_ iabPi CALL sARjxa,ZN,-Ts,ii7,zn ^ ' 

L''tf.7a ' . Y78Ya 

'lh80_ _ CALL SAP fX>jY'7',-RL‘*27V?) ‘ __ _ 
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c: 

ri 

c: 

c: 


SEO 

1 7P0 
7 in’ 
_ i 7?0 
[ \75d" 

1 7an 
i 7S0 
1 7 fed 
[_ '17'70' 

■ 1780 

n 'i 7 <70 
" 1800 
[ 1810 

' i 8?0 

f 18 30 

1880 
1850 
1880 
1870 
1880' 
18«>0' 
1 800 
1810' 
1 8?0 
F T 830 
"i 8U0' 
r ■■1 850" 
'■--1880^ 

“ "1870 

I 8 80" 

1 1 880^ 

“■'■?0O0' 
[■■■■?010 ■ 
? m ? oi" 
?0 30 
?0U0^ 
■?050 
?0h0 ■ 
?070 ■ 
?080 
I" P080 

■ "P100 ■ 

r" ■? 110' 

"■?i?0 

[ ?1 30" 

-'■?ia0 

I ?150- 

— ?lfc0" 
r' ?i 70' 
?180 
r— ?i80 ■ 
-■ P?00‘ 
[■■ ■PPI^" 
??? 0 “ 

12 

J ??U0 


CZSMP 

i t • • • 


,FTN 


1 , 


pb-FER-ai iptaaipa 

.i I.. I, 


PAGE 


T68‘ 


Z?«Z7 

call SAR {ZP.X?,PT,Z3,<3) _ 

CALL SAR (Y?,73,Hn,va,Zfl) 

CALL SAR (Z«,X3,GCLT,75,K5) 

'c all sap ' c va , as , ln, Vh, x6) 

"Aosi(<!iW)(h*¥*)»v#t*zs«zsAecs • 

« f A8* «0 + V8»y0*75«Z0*ECS)/AQ ■ 
Tt.SA (-HtJ-SOPT (Ba*Hf5-a,*C0/AO) ) 

■ At 1X8» r*)(0 2" _ _ 

■yEtYfc*T + V0 “ 

7E»Z5*t*70 

I T L H • A T A N ( 7 E A t C S / S Q R T f X E * * E ♦ V E * V E ) )■ 
■LNLP»nPW*4TAN?fye,XEl 

■'L iLHaDPP»LTLP _ 

CALL GT (LTLP,LNLP,1PN,U/V) 

'Ot (U-U0) /SF - '2 ~_ 

■ Vi (V»V0) /SF ■ 


I F C L » . r. T . ? V ■ W H I T E ( K B , 1 6 8 )" L TL P » L NL P 0 7 V 
FORMAT ( 2Fl?-,3,aFt?,A) - 


TT0' 

IT?- 

"C 


■K»l 

‘L) Tii^,'“ 

'00 175 ■i»r,8 

vj«i. 

'00 I70"jil7(i.' 
"VfFIa IJIaVJ ■ 

VjiVJiVT 

Ksn"*1 

UTfuTaU 


r 


r ■ 

L._. 


Tfl0 


:i 80 - 

C . 

"c ■" 

C 

C" ■ 


~TF "f I TT.EO.TVXM J«TK 
" I F ( 1 T T . e O'. 2 X ( 1 ) ■ ! 'f ;■ _* _ 

NsM*l^ 

M«0- ' 

no 180 1 i I7i P 

_XME An'C I ) aXME AN f i l VX ( I J“ 

Do 1P0 ‘ j«I ,TP 

'MfP + 1'” 

ACM) «A (Mj + X (T) AX t j) 

CONTIIYLIE“ 

XNaAI ■■ 

'aN'T^sW 

■RESOFa^xW ■ "" 

'■Ms0 ' 

"»fcsnF»xN-F ,0 ■ 

no 1 Q 0 liiVIP 

■ n 0 ■■ 1 8'0 j ■ 1 ; I p ~ " ■ 

M«ptL" 

A ( M l « A (P 7 - X M E A N ( I Ij X H E 4 N ( J ) / X N" 

REPLACE" XMfAN WITH" MEAN VECTOR, 
COMPUTE STANOARn DEVIATIONS ' 

' MA0 


A WITH covariance MAITRIX, AND 
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SEQ 

C7SMP 

«FTM Pb«FE8>ftl lPia4l2B PACE S 




00 200 I»1,1P 
XMFANdl «)(ME*Nf n/*N 




no ppp jti.XP 




M«M+1 . . • 


r PP90 


ArM)«A(M)/PESDF 


PTPI0 

POP 

CONTtMUe 


r psip 


DP PIP I«1.IP 




LTl.LFfl.I) 


F 


STOEV(n»OSO«T(A(LIin 



PIP 

cowriNUE 


1 PTSP c 

?lbP 

c 




c 

replace upper diagonal section of matrix WITH CORRELATION 

matrix 

?>I80 

c 



... 


Map 


papp 


00 P50 I»1,TF 


r poiu 


00 P5P J«l. IP 


papp 


MiM*1 


( pajp 


IF n-ji ppp.aap.ppp 


paap 

??0 

POSOaSTOEV f n*ST0EV{J) 


[ pasp 

' 

IF (POSO.tO.0.01 GO TO 230 


pafap 


A (Ml aA (Ml /POSO 


1 paTP 


60 10 PSP 


?a«0 

P30 

-AfMlap, 


r paap 


GO ro PS 0 


PSPP 

pap 

A(H)al, 


r psip 

P50 

CONTINUE 


P5PP 

C 



f P55P 


XNaXN-1 .0 


psap 

C 



L-?55P/ 

C 



PSbP 


KOEPal 


r psTp 


FINCa.5 


PS8P 


FOUTa.l 


r PSap 


tol»i -E-10 


PbPP 


HA»5TPaIP*P 


r pbip 


no PbP lat , IP 


pbpp 


cmap.p 


r ?b?p* 

PbP 

Continue 


Pbap 


LOOaLF taOEP. KOtP) 


r pfr ‘^0 


CfKOEPJal.P 


PbbP 


NVOam 


r PbT0 


00 ipp lal.IF 


PbPP 


IF (T-kDFPI ?70,3PP,P80 


r pb”® 

P70 

Lin«UF(I,KOEP) 


P7PP 


NVnaMVO*! 




GO 10 ?R0 


prpp 

pep 

LIOaLF (KOEP.M 


r ?t;<p 


NV0»NV0^1 


P7ap 

pqp 

FE INVO) aAfLIO) aaPaXN/n .•A(LlO)aaP) 


r PT5P 

300 

CONTINUE « 


?7bP 


OFaP.P L 

- 

f ?77P 


Lap 


?7ea 

C ' 

310 

C 

l.«L + l . 

call SUHROUTIne ¥n enter variable, calculate values to BE 

PRINTED^' 

P8PP 


CALL STEPRG (AnVFt.) 




- 



1 


! 


78 



2™ 


ORIGINAL PAGE E3 
OF POOR QUALITY 


APPLICATIONS DIVISION 


CZSMP 

^ t • • • 


.FTN ?b-PEB>St UiaaiSB PAGE fc 

I......... ..i ....I 1.. 

IF IFLAGI 320,fc?9,338 

EMTERiSMVE ^ 

7gO TP‘'!io0^”’ _ ~ 17 

■ FNitfiiKNTM _ ■ * II ' ■ . 

" ino«LF (KPtp'.FhEP) ‘ T7' II 

RFSI)F«XN-UF' ; _ 

' PtSSS* «N» (STriev (KbEP) ««2l* XtLD C>V ' ' JIHL _ 7 

"PtSMS*RfcSSS/RtSUF ■ -- - 

PS8U«SC1RT («E988'/XNtj 7" _ I ■" ^ ^ 7'‘ 

RFROFsOF" 

'■REGSStXN^fSTOtVfKtrCPrAA^f-BESSS' _ 7^ 7~II7I_ 7 

“RERMStHtGSS/REGOF 1 177 7 - - — - 

' fratici«rfgms/re9ms777 —1117 I ^_7_7_7777 11 ! 

" STEWR«SOP r (HESM8J * ' ’ 17 

lifMiJLT««osQHTrii-A(Lbonl_ 77171_77 __71 ^7 7177_1.-7.17 

InFiOF ’1 1_ 7_1_7_ _ 1 

■'i«bF«RtsnF _771 _ 7~1^ 1 .7 77 1- 17_ _ 1 

"IF (LR.Lt,?) GO TO '3T0 

'WRITE flP, 359V L.eNTFRVKAY.XMULTRVSTERR.PSSO ■ ~ - - 

"FORMAT ('I ' STEP NUMBER »,T3S;i3/T5, 'VARIABLE 'TAR, ' 0' ,T35 , 1 3/' 0 ' 

T multiple fl',T3S,F«.4/T5,'ST0, ERROR OF EST T31 , Fl 8,«/T5,'RF.810U 

■?AL"SAMPLE 'RMS' 0EV.',T31 ,F1?,4) 7‘ ' 

' WRITE' (LP. 3t>0) IDF ,HEGSS,ReGMS,FRAT10, IRDF , RESSS, RESMS 7 

" format' if 'W analysis OF VARIANCE '/TPS, 'OF SUM OF SOUARES ~ 

T ■ "mean square F-RAT10'/T13, 'REGRESSION', 3X,IaV3X, 7 

PFib 13 , SY , F ra'. 3 , 3X , f_l a . 3/T t 3, 'Rt SI DUAL' , 5X , 14 . 3X , Fl6. 3, 3X , F 1 « .Sill ' 

'"A'VARlAeLf IS'IN THE' EQUATION IF C(IJ IS LESS THAN OR EQUAL tO'0,0 


'■ 

[ 3?1 0 

" 3 ?? 0 ' 
[ 713 ? 30 
3?«0 
[ 3?S0" 

'3?fe0 
'■"■3? 7 0 

3?S0' 


NVI»0 _ _ _ 

■Nvo'e0"'' , "" ' 71 '1771 "7"^ ] 7 

'ACPHAiVMEAN'lKOEP)'" 

"oo"«h0 'i«iiip’' 7117' '71-17 '7'1__717 77 _ 7"77 

'lit»lf(T,'i 7 " 7' ;7_7_ ’ ^71 ! 

■■lF"(7iVOEPr'380,aA0,5<>0 _ __ : 

"LiD»LF(r,>oEP) _■ 7717_ 17 7! 1 1_ ’ 1- 

"GO'TO 400 '7 ] : 

LtnsLFfKOEP,!) 

7 Fl'i c ( ir, GT . 0 ; r ~Go "TO '» 1 0 ~ 2 11.7 1771.7 7. ' 

■'COMPUtX MULTIPLE REGRESSION' EOUAYTON"COEFFIC1ENTSV8 TO', error 
"ANO F'Tft REMOvErFOP' variables IN THE REGRESSION 

"n V I •nvTa i 

'pTNviV»STDEV'fK0EPr*'A{LinV7ST0EV('I) ' ” ^1 

n(NVn»(STERR/STOEVm)*OSORT(-A(Lin7XNJ 

FfNVIVVfBrNVTV/OINVm**?" 

■■■ALPMA»ALpHA-n(Nvn *XMEANn) _ 

'INENtNVlV«I ~ ■ ' 7 ] 

' GO ' TQ"43i7i ~ ^ 17_ , 1 7 

'a'VaPTAPLE IS'OUT 'OF^tHE'REGRESSION IF C;cn"JS'76REATElR THAN'JORI 

■'EQUAL TO 1 - - --- • - - - - 


'C0MPUTE"PARTIAL"C0»Rfl.*T10N COEFFICIENTS', TOLERANCE, 
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SEO 

3370 ) 
^33801 
j3'»(i' 
3UO'«' 

3«?^~ 
r 3'i3g)~ 
3uao) ■ 
3aS("' 

' '3U8P)' 
[ ■ 30 7 0 ' 

[" 3O<j0 ' 
" *?SMPI 

Cssin 

3S?{» 
[■' 3530 ' 
3S«f) 
[_" 3550' 
“ 35bO“ 
[""3570' 

'■■5590 


CZSMP.FTN 

I ) 

C F TO ENTER 

C 

ain 


?b-FEB>8i 

...1 1 ... 

FOR VANIaPLES out 


12I4AISB 

I. 

THE 


PAGE 


OF 




REGRESSION 


a?ipi 

'oivi" 

oa^ 

"iTsiT 
4 bo” 

c 

~c ~ ' 

a 70" 


NV0«NV0*l 

"NIEN'(NVb') St ^ [ H ~ 1 

ItbLtVtNvojsAcLiii I “ _,_lTr’lTZ 

RtNv6TsAfLtD)/DSORtfAfLli)*A(l,6Di') 

F E ( N V 0 ) s ( A f L 1 0 ) i * a * f RE S_0F . 1 , / ( A (iillA A ( 1,00 ) ■ 

IF (I -8 AY) abf),a3O,«60 . _ 

'if fCfin aao,<i«0,a50 

F'kA Y«F (NVir 

”6b ■ TO ftbo '_'7' ' ■ 

FKAYsFE (NVOl 
continue 


LA(l.lD)_**aL 


TP flR.LT,?) GO TO 310 
WRITE heading for COEFFICiENfS 


WRITE (LP,a70) ■■ 

format f/57X,TH;7?ri(72YHVARIA8rES'iN"M^^ 

IBLESNOt IN EQUAflbN/57X, lH7/6)(,«'MVARIAHLE,6)«,liMC6EFFictENT,2)(,10 

ERROR, a«, IIHF TO REMOVE . , 5V , flHV AR I ABLE , OX • 1 3HP ARTI AL CORH, » ’ 

‘^‘ TO ENTtR/57X,lM,) ” s'— 


?HSTO. 

'35*, RmtOLEwanCe, ax, J 0HF 


L 3500 

C 


: 

3800 
[_" '3810 

c 

PP1N7 THE REGRESSION ANALYSIS TABLE 


c 



3820 


WRITE (LP,«R01 alpha 


C 3<> J® 

as0 

FORMA t i 5 7 X , 1 H . / 1 0* , 9H iCONST ANT, 1 X , F 1 1 . 5 , 2H 

).24X,IH.l 

3800 


NOOS0 


[■ 3850 

R90 

IF (NVOJ 530,530,500 

. 

3880 

500 

IF (NVI.LE.01 GO TO 580 


1 38 7 0 


LNV*MIN0(NVI ,NVO) 


3880 

C 



r 3890 

C 

NVO AND NVI BOTH POSI T I VE, PRINT BOTH SIDES 

DF TABLE 

3700 

c 



i 3710 


00 520 Isl.LNV 


37?0 


WRITE fLP,510) INEN(n,B(n,0(n,F(I),NIENfn,Rtn,tOLrV(i).PEfiT 

[Z3I30 

510 

FORMAT (PX, 13, l*,F19.5,lX,Fll,5,lX,FtI,l,3H 

_,,7X,13,IX,F19,5,IX,F 

37U0 


113.8, lX,F12,a) 


r 3750 

520 

continue 


3T80 


NVIiNVI-tNV 


L 3770 


NVO»NVO-LNV 


3780 


NRO«l. NV 


[_ '37 90 


GO TO 490 


3800 

C 



L 

C 

Nvn ZERO, PRINT LEFT SIDE ONLY 


3820 

C 



r 3830 

530 

IF (NVI.LE.0) go to 590 


38U0 


00 550 I si, NVI 


[_ 38 50 


IIsi*NGO 


3880 


WRITE (LP,5401 INENHI), Bill) ,0(in , Fill) 


1 3 8 70 

540 

FORMAT f 8 X , I 3 , 1 * , M 9 , 5 , i X , F i 1 , 5 , 1 X , F I i . i' , 3 H 


3880 

550 

CONTINIJE \ 

- 

[ 3890 


NVT s NVT ♦ LNV 1 RESTORE TO NUMBER OF VARS 

IN REGRESSION 

3900 


GO TO 590 


( 3910 

C 



3920 

C 

NVI ZERO, PRINT RIGHT SIDE ONLY 
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?f.O 

CTSMP, 

• ..\ 

PTN Pfi.PEB.Bt 1PI44|?S PAGE A 





1 



.1 1 


5« iP 

c 






ShP 

PO SAP I»H,NVO 



1 

JOSP 


niI*K!GP 






w9nf (IP.S7P) NitNnn,R(in,TOLEv(n),FE(in 

. •• 

, 

i 

HQT;^ 

STP 

FORMAT (S71(,lM,,T>(,15,ll<,Ft9,5,l»,Ft3,«,tX,Fia,a) 





SSM 

COMTINUt 





S9P 

TF tl-MAVSTPl 31P,APP,f>P0 




<1 

f,PP 

aiwi tf. (LP.M P) 




API P 

h 1 P 

format f?3H SPECiFJEP STEP HEACMFD) 




ap?p 

6?P 

RMS«RSSn*iiEi.Taf n t) 



i 

an^p 


TF (TTT.FO.n MMTTE rLP»6S0) RMS 





^^P 

FORMAT ('PFAST-kEST FIT*,F7.t,» METERS RMS') 




Oft-rC* 


IF (JTT.FO,?) write (LP,6aP) RMS 





T>ap 

format t 'PMORT h-SOUTH FlT'.FT.l,' METERS RMS') 




a,’ TP 


WWITF (LP.A*;?) MVI , ALPHA, (INENd) ,0{I) , 1»1 ,NV1) 




ap8P 


WRITE MVI .ALPHA, {iNENm .P(I) , I«I ,NVI) 





bS0 

format nS.FlP.H) 




a 1 Pp 

fefcp 

CO^TINU^ 




an? 


riOSF ( lINITtJ, OtSPa 'PRINT') 




a 1 ?p 


STOP 'RESULTS IN CZSMP.lST' ' ' " 




ai HP 


ENO 




a 1 ap 


subroutine STEPRG {APVFL) ■■ 




a 1 SP 


LOGICAL AOVFl 




ai hP 


OOUHlE precision X,XN,A 




ai TP 


dimension trap), Aiea?), c(a0)i F(«b), FEtR0) 




ai j>p 


common * , HN, A.C.F.f E.L.NTGC, IP,NIP,FINC,FOUT,KAV, 

FLAG 

,KDEP,DF,TOL' 


a 1 9P 


LFflH, JH)tJ*.i*4NIP*l.((lP-l*al)*np«lx*?))/2 




a?pp 


VMlNtl,E?0 




a?i P 


VMAXI.J , 




a??p 


Kvlap 




as Jp 


KVD»0 




asap 


no HP Kt 1 , IP 




a?sp 


IF (K.EO.kDEP) GO TO 30 




a?fjp 


IF (C (A) ,r,T,0,) GO TO 10 




a? rp 


KVI*AVU1 




a?8P 


V»F (KV n ♦ (C (Kl *7, ) A«a 




a?op 


IF (V.CT.VMIN) GO TO 30 




aHPP 


VMTNbV 




aHi p 


k N ■ K 




a Hpp 


GO TO 3P 




apHP 

10 

A VOaHvO*! 




a Ha p 


IF (CIAI.EO.I.) GO TO 30 




aj5? 


V»FE (AVO)* (C (A).?,) aaR 




BHfcP 


TF (V.LE.VMAX) go to 30 




aTTp 


LKRil FfK.Kl 




'anap 


IF (A (LAK)-TOLT'10,20',20 




a 39?' 

- ?0- 

VH A E A V 




4t6P£^ 


K M A V aK 



i 

aaip 

30' ■ 

continue ■ * - ■ 




aap? 


IF'I.NOT.aOVFL) 'GO TO 100 " ‘ 




a a HP 


IF "(vmin-fouT) ap,a0,60' ' 




aaap 

dp 

IF fi-n sp.sp'.sp ■ 


' 


aasp 

SP — 

C fKMJ N1 IC (AMIni >9, 




aaibP 


KAVaAMIN ■ ■■■■ ----- - 




aar? 


FL AGa-1 , ■ ■ " - ... - .. . - 




aa»^P 


GO TO 8P ■ ’ 






• 
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2™!!! 


sto 

czsMP 

,FTN 




4490 

60 

TP (VMAX-PiNCl 

4SOO 

70 

CrKMAX)iC(KHAX 

[ as to 


KAVaKMAX 

aspo 


FLaG»l . 

I as JO 

RO 

CALL STEP 

asao 


DF.nF*FL*B 

L 


HE TURN 

asho 

90 

PLAGaO. 

[ asro 


RE TURN 

asBO 

100 

CONTINUE 

r as90 


IF (OF.GT.l.l 

afoo 


GO TO 90 

r afcii? 


ENO 

4hP0 


subroutine lgr 

1 atijo 

c 


Ob Hit 



i afeso 

c 



Obbdt 

c 

FNVALl « 

f 4h70 

c 

1 

4h«*0 

c 

1 

[ 'Ih9ig 

c -■ 

i 

ayoo 

■ ■"C 

t 

9 

r aTio 

C il 

a a a a ^ . a ^ a . a . . . . 

a 7;?0 

c 


r a;j0 

c 

VALt - 


i...; 


12I44IP6 


KRESLT 


PNVAL2 - 


_FNVAL5 - _ 

I FNV AL4 . 

„.l J. 


• 4 » - •— ii,-r U * 


r_'aTB(* 

insvi 


j. ijit^ ?. V-* L . • y * w 3 .. . • . 


G1V|n>0uS' vaTU6S ■T'VAlF 4 or APRAT_yj AND _ 

^ ?. .“p./i-OR A“**TL.Fy.j.. that_ 


" aiRO 

C 

4N0 GIVEN A Value xintrp i xi ) that is adjacent to 

VALl •> 4 

r 4 7 90 

C 

then this ROUTINE RETURNS A VALUE, XRESLT# 

THAT 

IS 

THE LAGRANGE 

aRoo 

C 

INTERPOLATION OF XINTRP ABOUT FNVALl •> 4. 





r a«i 0 

c 






aspo 

c 






1 as JO 

c 



* 



aRoo 


real v( 41 , Fvrai , xreslt , xi, part(4) 





( a 8 S ;1 


real a.h.c.d 





4860 

c 






r 4 8 70 


A a xl • V(l) 





a 8 80 


B a XI - V(?) 





4890 


C a XI » V(31 





4 900 


0 a XI <• Via) 





[ 49J0 

0 






49?0 

0 

TyPt 500,A,B,C,0 





! a 9 JO 

0 

500 F0»HaT(4F1?,6 





4940 

0 

TYPE 500, V 





r 4 9S0 

D 

TYPE 500, FV 





4960 

n 

type 500, xI 





[ a970 

c 






4900 


PART(l) a f B a C a D ) 





[ 4 990 


a / (( vm • V(2)) a ( V(1J - V(3)) 

* 

( 

vm 

- Vf4)) ) 

SOOO 


♦ a FV fi I ' 




- 

r SO10 

c 






SOPO 


PARTf?) a ( A a C a 0 ) 





[ sojo 


a / n v(2i - vnn * ( v(?) • vtsn 

* 

t 

V(3) 

- vian ) 

SO40 


♦ a FV f?) 





... 


- 
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SEQ 

CZSMP 

.PTN 26-PEB»0t 12144128 PAGE 10 



• • 1 • 


C 

P4RTf3) • f 4 * B * 0 ) 




1 SPiro 


* / rf vf3i - vnn * t vfS) - vfzn * f vtsi - 

V(4)J. 

J 


S080 


♦ * FV (3) . • ‘ ■ 




50 90 

c 

s’ ‘ 




5100 


9*97(4) ■ r A * B * C ) 




f 5110 


♦ / (( V(4) - vein * (' Vtfl) . Vf?)) 9 ( V(4) - 

vesn 

) 


51 ?0 


♦ * FV (4) 




L 51 J0 


X(iE8L,t ■ P49T(i) <* P4Rf(2) ♦ PARTti) ♦ PART(4) 




51«0 

c 





□ 5150 

D 

TYPE 500.P4RT 




5150 

0 

TYPE 500.XRE.5UT 




r 51 70 

C 





5180 

C 





1 5190 


RETURN 




5?00 


ENO 




! 5?10 


SUBROUTINE STEP 




“ 5??0 


OnuULE PRECISION X,XN,A 




f"5?30' 


OIMENSION Y(40), 4(840), C(40), P(40), PE(40) 


DF7 


5?U0 


COMMON" y.xN, A, C,F,FE',l.,NT6C, IP, NIP, F INC, FOOT, KAY, FrAiS, 

KOEP, 

TOL 

[■■■ 5?50- 


"LFfir,"Jx)«JYiiy*NIP*l.((TP-lX9l)*(IP-IX92))/2 




5?80' 


L8K«"tF(K4Y,K*Y) 




\ 5? 7 0 


M»0 


i 


■■■■5280 


no 70 i»i,iP 




f- 5?90 


00 70 J«t,IP 




■ 5300' 


MiM*l 




(■ 5310 


-"IF- '(I-KA-YJ- 10-,-r0V20“ 




5330 

10 

■ LIKiLFn.KAY) 




r 53 30 


'GG TO 30 




■ 5 3 a 0 

?0— 

LTK«LF(KAY, I) 




n‘5350 

•^•"^10 

“IF CJ-KAY) 40,70,50 




~"‘5 380' 

" 00 

LKJ«LF|.I,K*Y1 




r '5 3 70 
"5380 


GO TO 60 , 




50 

LKJ«LF(K4Y, J) 




r" 5390 

• "60 

" 4 (M) «A (M) iA (C IKT*A fll<'J)7A(^ 




5000 

- '7 0 

CONTINUE 




f"' 50 10 


"00 110 it'i.iP 




5't?0 


IF (I-K4Y) 80,110,90 




Q” 5030 
■■ 5000 

■•■"50 

L1X»LF(I ,KAY) 





GO TO 100 




f""5O50 

"“'90" 

LIK«LF(8AY,n 




" '5080 

100 

A(LIK)«4(LlK)*FI.AGyA(LKK) 




{ 5070 

" ria"" 

"■ CONTINUE 




“5O80 


.00/AaKK) 




j 50 90 


RETURN 




"■5500 


end 




I 


1 : - 


1 1 


1 


1 


1 . 
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CZCS NEAREST NEIGHBOR RESAMPLING PROGRAM 



sen 

C7SNR, 

FTN/5 2O«0E8«SI 12146108 PASS 1 










101 

C 

C7C3 NEAHtST RE8AHPLER 




2d 






501 

C 


- 




C 

AOAPTeO BY GLEnm OAVIS ON SEPTEHBER 21.1480 PROM .LN0 




^01 

C 





fcOI 

c 

' 




701 

c 





SOI 

c 





9M 


PARAMEtER NS(.«34«N6La620.MAXOLV*i6.eyTRYTB0.URDW>b«2 




101? 


iNTEGtRAP WC.wR 




1 IP 


lNTERErt*2 I0FE{NSL),10(5), CVAU, START, FIRST 




1 ?0 


LOGICAL*! LHP(?0481 , JttUFfNBL.NSU ,ST (10) ,STO(10i ,HFN(18) 



t 50 


LOGIC*L*l PRNTJH, PRNTHD 




140 


LOGICAL*! DA(101,TM(I0) 




150 


lMTtr,ER*2 CRN, IV (MAXOLY) .NftRMVS.MVTYPE.lNSklP.OUTSKP 




1R0 

c 

NfiLTS IS NUMfiE.7 OF 8LOCKS TO SKIP ON SHORT READ 




1 70 


OtMtNSIOM DHVAL (6) .hCVAL (6) , OHfS), bC(6) 




1 80 


REAl. UTSC.LNSC 




1 90 

C‘ 





?00 

......p- 

Call ASSIGN! 4, »CTCS,L8T») 




210 

c 

X 




2?0 


CALL liATEtbA) 

; 



?50 


CAIL TIMp(TM) 





c 

' 




250 


WRITF (5,201) OA.TM 




200 

“O"' 

WBIIE (4,20) OA.TH 

, 



270 

20 

FORMAT (*0C7CS NEAREST RESAHPLER V 1..0 * , 4X . 1 0 A 1 . 4X . 1 0A 1 /) 



280 


NSLMtNSL»l 




P90 


NwLsNhL/2 




500 


N8TH«NeL *512 




510 

c 





320 

c 





5 50 


iall»«ll 




340 


WK11F (5,30) 




550 

30 

FORMAT ('SINPUT IMAGE TITLE t ') 




500 


REAO (5,40) ST 




570 

40 

FORMAT (10A1) 




580 


CALL GMEN (ST.HFN) 




390 


Call ASSIGN (2,HFN) 




400 


read (2,50) 10 




410 

■50 ■ 

format MX,5A21 




4?0 


read (2,801“ Nlf7Ntt',N0P, IIP, NlPPVlEb.ILb' “ 




4 30 

00 

format (715) ■ 




4 4 0 


read (2,70) OELPVOEUC, UTSC.LNSC 




450 

T0- 

FORMAT (4E12.6) 




400^ 


CLOSE ( UNIT * 2) 




4 70 


CALL PI (lO.IPNl 




480 


NHIsNIt*N8P 




'4 90' 


NEt.® NBL/NBP 




500 


TF (MR.NE ,21 1G0“ Td TT« 




510 

77 

WRITF (5,801 ' 

- 



format ('SWHTCH DHIVC ~7 *) 


SSO PfAO f5,R0) noR 

SaPI ■ RP> FORMAT (TI‘5) 

' SSP''~C 

■ 5»)PI WR I TE 'rs, IPfl) 
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SEQ 

CTSNR, 

FTN13 96-FEB«8t 12146108 

PAGE 

2 



1 





571! 

100 

FORMAT f*#nuTPU7 IMAGE tlTLE 7 *) 




S 8 CT 


9EA0 15,401 8T0 




S9PI 


CALL GHFN fSTO.HFNl 


' 


hPiUI 


CALL ASSIGN (2,MPN) 




6t 0 


REAO (2,40) 10 






READ (3,601 NE,NL,N 8 P,IP.NPP 




630 


READ (2,701 OX.OV.FLT.FLN 




6«0 


Nl A bNSL/3 




6?(0 


CLOSE ( UNIT • 21 




660 

C 





670 


IF (HO0(NBP,2) ,EQ.O) (iO TO llO 




6B0 


INSKIO ■ 1 




690 


QUTSNP a 1 




700 


NPKMVS a NBP 




710 


MVTtPe • RYTBYT 




7?0 


GO TO 130 




7 30 

c 





7U0 

110 

IN 8 KTP • 2 




7S0 


OUTSKP a 3 




760 


NHRMVS a NHP/3 




7?0 


MVTVPE « WROWRO 




■'7 80 

C 





790 

120 

CONTINUE 


; 


B00 

C 





610 


WRITE(5,801 




M?0 


RF.AO (5,901 lOOR 



! 

830 


nho«ne*nbp 




flU0 


NBTMi 1034 



F” 

8S0 

c 



• 


. 660 


MONaOV/3, 




670 

' 

call GT (FLT,FLN,IP,XF0,VF0) 




660 


CALL G7 (LTSC,LNSC, 1P,XSC,VSC) 




890 

c 





900 


write (S, 1701 




910 

1 70 

FORNAT {'SBEGINNING ROW AND COLUMN 7 ' 

) 



930 


PF4U (5,1801 IHH.IbC 




9 50 

180 

FORMAT (31101 




900 


IF (lHH.Nt.01 GO TO 200 




9S0 


IHBaj 




960 


IHC »1 




970 


IEP»NL 




980 


IFCCbmE 




990 


WRITE (5,190) I8R,IBC.IER,1ECC 




1 -Tiwpi 

0 

WRITE (4,1901 T 8 R, IBC, UR, lECC 




1010 

190 

FORMAT ('0OEFAULT VALUES! '.aXS/J 




1 


GO TO 330 




1 0 50 

300 

WRITE (5,210) 




100 0 

PIO 

format ('SENDING HOW AND COLUMN ? •) 




10S0 


REAO (5,180) IER,IECC 




1 060 


IF (lER.GT.NL) lERaNU 




1070 


IF (IECC.GT.NE1 IECONE 




1 080 


WRITE (5,60'1 ' 


- 


1090 

D 

WRITE (4,601 




1 1 00 

??0 

NOE a 100 




1110 


IRMNf ILO+l 




1 1 30 


TRMXa ILO^NIL 
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C7SNR 

,^TN13 06-FEB-M t?l46l0B PACE 3 








.. J 

t . 



1 1 SCI 


TCMN» Jt04l 





1 1 40 


icMx. no*Nte 





1 1 ^*0 


CALL ALTPHU,40) 





1 t bVI 

C 






1 1 T0 

C 

V ■' 





1 1 «« 


CALL tMOPtN (t,8T,,lI0R,NHI,NlL,»P*) 





t 190 


CALL IMnPFW f0,STO, , 100K,N80,NL,'M») 





1 ?0ii 

0 

CALL ASSIGN ( 3, 'CZCS3.LST*) 





1 ?1 0 

C 






1 ??0 

c 






t?J0 

c 

"WRltP 19, 1941 XAO, rfO,XSC,YSC 





i ?40 

C14« 

FORMATI * XFO VFO X8C 

Y3C 




1 ?sn 

C 

♦ ,/,9F15.4 1 





igbO! 

c 

WRTTF (4, 148) tNSKlP,OUTSKP,NBRMVS,MVTVPE 





1 ?70 

C14H 

FOPMAtf ' INSKIP OIJTSKP NflRMVS HVTVPC »,/,9T® 1 





1 ?fl0 

C 






1 ?90 

C 






i 300 

?30 

1EC»I8C+N0f-1 





t 310 


IF' (ItC.Lfc.lECCr CO TO 240 ■ ■ ' 





1 3?0 


iEC»iecc ^ - - - - 





i 330 


■ NOF*lfcC-lHCYr 





1 300- 

■ ■?40 

NOL*Ifc»-IBB*l ■ , ■■ " 





1 350 


IKCPbUC*! 





t 3b0 


K0» ( TBC-1 ) *NflP 





1 370 


N0»K0/51? 





1 380 


TFfNO,lT,0) NO»0 





1390 


KO«ko-NO*51? " ■* 





1 400 

c 






1910 


THESE values are swath tlHlTS 





19?0 

■'C 






"19 30 


OL* ■' XFrt-KSC 





1 940 


DLT • rpo-YSC 





1 450 


XMINbOLX* ( IBC-,51 »0X 





1 Ot»0 


>MA**xMIN*fNOF:-n«OX 





1470 


YMAX rOLT- ( I0P-,51 *DV 





1 4R0 


Yf^IN*YMAx - (MflL-1 1 BOY 





1 490 

. c 






1 S00 

c 

THt 21 TfcRM POLYNOMIALS DEFINE BOTH U AND V ADJUSTMENTS, 




1510 


SINCE THIS SCHEME WORKS ALONG A SINGLE OUTPUT LINE, 

THE 




15?0 

■'■' C 

VERTICAL COMPONENT NEEDS TO BE CHANGEO ONLY WHEN A 

NEW 




1530 

c 

OUTPUT LINE IS CHOSEN, P2 EXTRACTS FROM THE POLYNOMIAL 




1540 

- C' ■ 

COEFFICIENTS THE HORIZONTAL ADJUSTMENTS, OR SO THE 

STORY 




1 550 

■ C 

GOES, " ■ 





I5b0 

r 






157 0 


CALL P2 fXMIN,XMAX,yHAX,y,DX,DY,OY! ) 





1580 

c 






'1590 

c 

WRTfE ■C9, r90r YMAX,xM1N,YHAx',YMIN ' ' ' 





1600 

"Clff0 

FOflMATr'' "XMAX ■ XMIN “ ‘YMAX 

“ ■ “YMIN 




1610 

--C' ' 

♦ ■ 77,4Fr5,9') ■ ' 





1 600 


WRITE (4,142) OX,nY,nLX,OLY ■ ' - 





1630 

CH? 

■FORMAT I ' OX - { OY ■' DLX '■ 

♦ ,779F15,9 ) V ” 

OLY 

#" 



1640 


■ 




1 650 

c ■ - 

WRITE (4, )96V TH,Y,0VI 





1 660 

C146 

FORMAT.f * TH ' ' ■ Y DYt 

■ ■ ■ ■ ‘ t 




1670 

c 

♦ ./.9F15,a 1 ■ 





1680 

C ‘ 
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SEO 

CZ3NP. 

PTN73 P8-PEB-81 12196108 PAGE 

a 



i 



. 1 . 

1690 


09NIIB9 



1 T00 

0 

WRITE Ca,60)IEC 



ITIP 

c 




t 7?0 

c 




1710 

c 




1700 

n 

KH*nO ■ 6 



1 7?)fll 

0 

KMIN ■ 9999 



I7fc0 

D 

KTIMKW ■ 0 



1 770 


00 600 J9l ,N0t 



1780 

c 




1790 

c 

PI calculates the 0TH THROUGH 5TH DERIVATIVES 

OF BOTH THE 


1 800 

c 

ROW ANO THE COLUMN TERMS WHICH GIVE THE CHANGE FROM ONE 


1810 

c 

EVALUAITON op the polynomial to the NEXT. 8V 

ADDING THE 


1 8?0 

c 

CHANGE TO THE PREVIOUS VALUE OP THE FUNCTION, 

A LOT OF 


1810 

c 

calculations ARP AVOIOEO. THE 0TH DERIVATIVE 

IS THE 


1 890 

c 

FUNCTION ITSELF. THE DERIVATIVES ARE WITH RESPECT TO THE 


1 8<j0 

c 

EASTING OR ALONG ROW DIMENSION STARTING FROM 

THE RIGHT SIDE 


1880 

c 

OF THE SWATH. 



1870 

c 




1 880 


Call pi (V,RL,CL,0R,0C) 



1 890 


OClila 0C(1)-1 



1 900 


wc« octn 



1910 


W»* DR(1) 

; 


1 9^0 


ICLa CL 



1930 


IHL» RL 



1990 


IF tJ.Eo.n ilnbiri. 



1 9S0 


NPN9WC-ICLA1 



I960 


NLNbWR«I»L*1 



1970 


IF (NFN.GT.NELl STOP 'INSUFFICIENT ELEMENT STORAGE' 


1980 


IF (NLN.GT.NSL1 STOP 'INSUFFICIENT LINE STORAGE' 


1 990 

' c 




?000 

c 




?010 

c 

INVOKE HLI J8UF 



?0?0 

c 

ASSIGN 99901 TO 199900 



?0 10 

c 

GO TO 99900 

* 


?090 

C99901 CONTINUE 



?0S0 

C*«**i 




?080. 

c 

PROCenURp TO FILL JHUF 



?0/0 





?080 

c 




?090 

C99900 CONTINUE 



?100 

as0 

IF (IHU+NSLM.tE.ILN) GO TO 280 



?1 10 

c 




?1?0 


VVAL 0 Y ♦ f ILN - IRL 1 * DVl 



?1 10 


CALL P3f VVAL, RLVAL, CLVAL, ORVAL, OCVAL ) 



?1 90 


CVAL » OCVAL(l) +10 



?150 

c 

ICl VAL » CLVAL 



?lb0 

c 




?\ 70 

c 

ICLVAL IS PIXEL CONTAINING LEFT EDGE OF INPUT 

WINDOW. 


?180 

c 

SU8TR4CT ONE TO SKIP OVER ONLY THOSE THAT ARE 

TO THE 


?1 90 

c 

LEFT OF THF EDGE, SUBTRACT ANOTHER AS A MARGIN OF SAFETY. 


??00 

c 


- 


??10 

c 

ICLVAL » ICLVAL • 2 



???0 

c 




??10 


LC«H0D(ILN-1 ,NSLJ»1 



??90 

c 
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seo 

CZSNR, 

FTNI3 P6-FEB-81 l?ia8i0B PACE 5 



1 


f 


P55 

IF(LC.GT,0) GO TO 256 


??feB 


LC» LC+NSL 


??7PI 


GO TO 255 


??PP 

?5b 

CONTINUE , . ■ * ■ 


??9?l 

C 

V * * ' 




STAMT « MAX( 0, CVAL « lEO - ( NBL/NBP 51 




NHlTS* START »NMP/512 


?^?PI 


IF (NBLTS.lt. P) NHLTS*P 


? ♦SP 


IRMILN-TLO 




IF (IBN.Gfc.l .AND.IBN.LE.NIL) GO TO 26B 




CALL HOWE (0,0,l.HF,?,102A,?) 


? 


(10 TO 270 


??70 

Pb0 

CONTINUE 




CALL IMSWRT (1 ,N8LTS,NHTR) 


?^<7PI 


CALL IMRt'AO n,I«N,LeF5 


?aw'^ 

?70 

NBVTSA NBLTSA512 


?«tP 


Ki 5TABT *NBP-NHYTSaI 



?7S 

CALL MOVE (LaE(K5,l,J8UF(l,LC5,l»E'BL,»5 



C 



?(IUCI| 

.... -J. 

The input SPACE IS SKHiEO IN RELATION TO THE OUTPUT 'SPACt;"' ' 


■ ?assi 

C ■ 

TMf OEFSETS calculated here allow the INPUT SPACE TO BE STOREO' 


?<lfeC1 

c 

AS A rectangular AUtA. WHEN IT IS INOEXEO INTO. THE OFFSETS 


?fl70 

~c 

ARE iiSEO TO SELECT THE PROPER AREA, T 


?««pi 

c 



P«90 


I0FF(lC)« STARTaNBP" " 


PS00 


iln#iln*i 



c 

IF (J.EQ.H KHIN « MIN( START, KHIN 5 


?5?0 


GO TO 250 


?5J0 

c 



?«!a0 

P80 

CONTINUE 


?5S0 

C 



PS60 

C 

GO TO 99901 


?S70 

c 




C99<>Pil 

GO TO 199900,(99903) 


?S'50 

c * • * * * 







Pbl0 

. c 





IF (Nn«512>N8TM,GT,NBO)NBTM»N80»NOA5l2 


' ?630 

'C ■ " 



?ba0 

c 





CALL IMSHHT (2,N0,NHTM) 


?bb0 


CALL IMRfAO (2,0RN,LBF) 


?h70 

- (. - 



?bH0 

-15 

IF ( HOD( .1«1 ,501 ,LE,0) WRITE(«, 40985 


?b90 

DA098 

FORMaTT'I *', ' - "LINE - HC "Id - 'NPN~ " . WR IRL», 


?700 

0 ♦ 

- A NLN STARTS » LINE* 


■ ?710 

U ♦ 

,TS0, ■ » 1 k WC LC IOFF(LC) WR»5 


?7?0 

-0 

W w T f E ( 4^ 4 0 9 9 )- J , W C ,"W R , I C L; I » L V N P N , N L N , S T A R T , L C " 


■ P7 30 

- OU0«»q- 

FORMAT (RIB ") ' ' . - - . - 


?7fl0 


KOOT« K0+(N(iE-T5*NBP+t 


?7S0 

■ D - 

IF -(M00(KT1HER, MSL5.lt, 15 PRNTJB ■ .TRUE, 


■?7h0 

... 

KTTMER ■ KTIMER ♦ 1 " ' . 


?770 

- c - 



?760 


no 400 I»1,M0E 


” ?790 


R»" DRd ) - “ 


?PP0 


WRi R 
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1 1 • i • 





CALL HOVE (0,0, IVjP.NBP, 1) 


PB?«t 


lEfWR.GT.IRMx.OH.WR.LT.lRMN) GO TO 3RR 



0 

tF(WR,6fc.tHHN.4ND.ttR,LE,lRM*)G0 TO 382 



0 

W»ITE(fl,5805 I , WR, IHMN, IRHX 



038PI 

FPBHATt • ROW VALUE OUT OF BOUNDS I I , WR , XRMN, JRMX » , 4 16 ) 



0 

RC TO iRR 


?S7(? 

C 



?8fl7l 

382 

LC«MOi)fwR-l ,NSL)Al 


?890 

C 



PQHKf 


C" OCM) 


?0|fl 


WCb C 


pqpol 

C 



P'J?!?* 

C 



?9ttP) 


TF (wr .LT. lCMN;OB.WC.r,T,ICMX)GO TO 3RR 



0 

lF(WC.GE.ICMN.AND.WC.Lt,ICMX)GO TO 384 


PPbH 

0 

WHIie(tt,5831 I.WR,ICHN,1CHX 


?9?a 

0385 

FORMATf * COLUMN VALUE OUT OF BOUNDS 1 I , WC , ICMN, ICMX ',4X6 X 


?9«01 

0 

GO TO 399 


?qq(|i 

C 




384 

K« WCbNHP.IOFF aCX At 


3»1 (? 

C 



3P?PI 

D375 

KOELT • KPMEV - K 



C 

WRITE (4,4(395) 1 ,K, vie, LC,idFF(LC),WR, KOELT 


i'Aua 

C4(i19S 

FORMAT(TH0,7I7 ) 



0385 

KPRtV » K 


3'’i<S(’l 

C 



3?i73 

0 

IF (K.LF..N8L.AN0.K.GE.I ) GO TO 388 


3?I8(" 

r. 

KtfiRCT • KERRCT ♦ X 



C 

IF ( M00( KERRCT, 100) ,GT. 10) GO TO 599 


III?!!* 

0 

WRITE (4, 586) J,1,K,WR,WC,I0FF(LC),LC 


IttP 

' 0386 

FORMAT!' K OUT OF BOUNOSI OUTROW OUTCOL ' 


31?'^ 

0 

♦,'K WR WC TOFF(LC) LC' 


31 3(?i 

0 

♦ / 16X,7(I6,?X)) 


3 ) aPi 

n 

GO TO 399 


313CT 

c 



31 fePi 


IF(K.GI .Nbl 1G0 TO 399 


31 TPI 


IFfK.LT.noO TO 399 


31 81^ 

388 

CALL move CJBUF (K,LC) , INSKXP, X V , OUT 5KP , NHRMVS , MVT YPE) 


31 9P 

• 

k»k»nhp 



C 



3?1<3 

59<» 

call MUVEdV, INSKXP, LHF(KOUT) .DUTSKP.NBRMVS, MVT YPE) 


3?8i8 


KOUT* KOIJT-N0H 


3?3C» 

C 



3?«C5 


no 402 J0«1 ,5 


325(3 


jp» jn +1 


328C3 


0C(J0)B DC (JO)+OC( JP) 


327P 

402 

0R(J0)8 OR(JO)aOR(JP) 


328t3 

C 



32953 

D 

IF (.N0T.PRNTJ9) GO TO 400 


35(3(3 

0 

prntjh 8 .false. 


331(3 

D 

PRNTHI) 8 .TRUE. ( 


352(3 

0 

KMIIV 8 ( lOFF(LC) / nAP ) - 5 


3 3 3(3 

D 

JPLUS 8 LC ♦ NSL •! - NLN 


33(H3 

C 



335(3 

0 

no 4083 KINDX 8 LC » NLN ♦ 1 , JPLUS 


556(3 

0 

KWHC 8 MOUfKiNO* - 1,NSL) ♦ 1 




- 


I 
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0 

MEMO • ( lOFPfKWHC] /M8P • NNIN ) * 4 6 1 


■578?) 

0 

KIND ■ ntXlKENO.n 


r 35951 

D 

KFMO ■ MfNf KENO,P5) 



n 

IF (.NOT.PBNTHO) t.O TO 4083 • ' ’ 

, 

3U> ft 

D 

PRNifHP ■ .FALSE, 


Sapfi 

0 

IMP* ■ 1067 fKWHCI /N8P 


5a in 

c 

INPX ■ IP6F (XWHCl/NbP ♦ j ‘ 


5aan 

0 

FIRST ■ K /NPR « KBANO ■ 1 


5850 

0 

w»ne(3,408n irl.inox, first 


308(? 

O80B1 

FOHMATL'I INPUT 8“ACE STARTS AT ROW *,14,* COLUMN *,I4, 


5870 

0 ■ ■* 

/• 0I5PLAVE0 FROM COLUMN »,I4) 


5880 

0 

WHITE (3,?01 OA.TM 


3890 

c 

IF (WENn.GT.flB.OR.KEND.Lt.P) TVPF *, » KENO ■ »,KEN0 


.5S00 

c ♦ 

lOFP (LC1 ■ ' , lOFF fl C) , • IC *,LC,» KMIN *,kMlN 


5510 

08083 

WRITE (3,8085) (JOUFriNOXiKWHC), INOX* FIRST, 156 * FIRST, NBP) 


35?0 

08085 

FriHMAtf T <KENO >,3008 ) 


3550 

0 

WR?TEf3,8088) IBC 


3580 

08084 

F0RMAT(»1«,* OUTPUT SPACE*/ * LINE PIXEL* / 5X,la) 


3550 

■ c ■ 



" 3580 

800 

CONTINOfc' 


■■ ■ 35 70 

■ ■■ r. 



3580 

... 

k£NO ■ NOE « NBP ♦ NBP 


■■ 3590 

0 

KENO • MIN( 180, KEND I 

; 

■ 3600 

0 

KENO ■ ( KOUT ♦ KB*NO »i ♦ KEND ) ■ NBP 


3510 

' 0 

WRITE('3,80H0i OHN, (LBF(IN0X1, IN0X8K0UT ♦ 5 ♦ NBP, KENO , NBP j 


3680 

08080 

FORMAT ( * *,18, 3008 ) 


■“ '3630 

“c 



~ '3680 


CALL I**WRIT' (a,ORN,LBF) 


■ ■ 3650 


WRITE (5,8901 OHN,18C,lEC,OR(n,OC(n 


3660 

'8'90 

format ('♦*,3I5,?F10,U 


3670 


OHNiORN+l 


3680 

500 

Y«V*OVl 


3690 

■" C 



3700 


IF (lEC.EO.IECCI GO TO 510 


'3710' 


IHC«IEC*1 


■ 37?0 


no TO ?30 


" 3730 

. c 



“ 3780 

" 510"" 

CONTINUE 


■■■ 3 750 

0" • -- 

CLOSE T uNrr«8, nrsp '«'*PRiNT* T' "" 


“3760 

•" p- 

CLOSE 1'UN1TK3,01SPo*PRINT* I '■ - - - 


3770 


STOP 


3 7 80 

C 



3790 

C 



3800 


End’ ■■ ■■ 
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